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Within QED, we examine several issues related to constructing a parton-model-based QCD 
transport theory. We rewrite the QED analog of the parton model, the Weizsacker-Williams Ap- 
QQ ' proximation, entirely in terms of phase-space quantities and we study the phase-space photon and 

^\ . electron densities created by a classical point charge. We find that the densities take a distinctive 

0^ ' "source-propagator" form. This form does not arise in a conventional derivation of the semiclassical 

transport equations because of the overuse of the gradient approximation. We do not apply the gradi- 
. ent approximation and so derive the phase-space analog of the Generalized Fluctuation-Dissipation 

^ ' Theorem. Together, this theorem and the expression for the phase-space particle self-energies give a 

, set of coupled phase-space evolution equations. We illustrate how these evolution equations can be 

used perturbatively or to derive semiclassical transport equations. Our work relies on phase-space 
propagators and sources, so we describe them in detail when calculating the photon and electron 
, phase-space densities. We use these tools to discuss the shape of a nucleon's parton cloud. 

• PACS numbers: 24.10.Cn, 25.75.-q, 12.38. Mh 

^ ' 

o 

O . I. INTRODUCTION 

oo 

Primary hadronic collisions in a typical nuclear reaction at RHIC will occur at y/s ^ 200A GeV. Such a collision is 
so violent that the partons, i.e. the quarks and gluons comprising the hadrons, will become deconfined. With hadronic 
densities exceeding the inverse volume of a typical hadron, the partons will remain deconfined and are expected to form 
a quark-gluon plasma (QGP) [0-^. Since transport theory descriptions of nuclear collisions have proven successful 
^ at lower energies, it is natural to attempt to describe the time evolution of the QGP using a transport model derived 
^ from QCD. A transport model would describe the time evolution of the parton phase-space]^ densities throughout the 
collision. The procedures for deriving semiclassical transport equations using time-ordered nonequilibrium methods 
are well developed I p rp l - In fact, there have been several attempts at constructing a QCD transport model based on 
these procedures |10p|l2[|, but each of them have their problems. Chief among these problems is that one either treats 
the soft long-range phenomena (in the case of or one treats the hard short-distance phenomena (in the case 
of [^0|), but never both in the same framework. Normally when one discusses transport, one assumes a separation 
between the interaction and the kinetic length scales. If one relaxes this assumption then one may be able to treat both 
hard and soft modes on equal footing. We have not done this for QCD, but we have made several steps toward doing 
the analogous thing in QED. Our techniques also allow for a simple connection with the parton model. This paper 
consists of three parts, each one using QED to describe different aspects of the problem of constructing a partonic 
transport theory. In the end, we use our accumulated insight to discuss the shape of the parton cloud of a nucleon. 

Before outlining the paper, we must say a few words about our formalism. In the first two sections, we use 
Feynman's formulation of perturbation theory. In Feynman perturbation theory, one specifies the initial and final 
states of a reaction and calculates the probability of going from the initial to final state. Thus, it is the appropriate 
tool for calculating observables for simple processes (such as exclusive cross sections) . For this reason, we use Feynman 
perturbation theory to illustrate how the phase-space sources and propagators work and to calculate the reaction 
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probability for some simple processes. We show that, in Feynman perturbation theory, the particle phase-space 
densities have a "source-propagator" form. Namely, the densities are a convolution of the probability^ to create 
a particle (the source) with the probability to propagate from the creation point to the observatio n p oint (the 
propagator). We are not the first to consider writing transition probabilities in phase-space: Remler discusses 
simulating many-particle systems in phase-space. Remler's work is not immediately applicable to partons because 
it only applies to particles with large mass. For more complicated processes, i.e. when we only know the initial 
conditions, we must resort to time-ordered nonequilibrium methods. In the last two sections we use time-ordered 
methods to derive the phase-space evolution equations, to derive the Generalized Fluctuation-Dissipation Theorem 
and to discuss the parton distributions of a nucleon. The Generalized Fluctuation-Dissipation Theorem shows that, in 
time-ordered field theory, the particle densities also have a "source-propagator" form. It should come as no surprise 
that we find similar forms for the particle phase-space densities since both formalisms are are equivalent descriptions 
of elementary processes. For simple tree-type processes in the energy-momentum representation, the moduli of 
the Feynman and retarded (or advanced) single particle propagators are the same, so one can rewrite the reaction 
probability in terms of either In fact, both sets of Feynman rules are special cases of the contour Feynman rules 
in Appendix ^ 

In the parton model, a cross section is a folding of the Parton Distribution Function (PDF) with the cross section for 
the partonic subprocess. The QED analog of the parton model is the Weizsacker- Williams Approximation since 
a cross section in the Weizsacker- Williams Approximation is a folding the effective photon distribution with the cross 
section for the photon absorption subprocess |14-16|. In Section ||, we write the Weizsacker- Williams Approximation 
in phase-space in several steps. First, we write the reaction rate density for our "partonic subprocess," namely the 
reaction rate for absorbing a free photon. By writing this rate in phase-space, we also illustrate how we do our 
momentum-space to phase-space conversions. Next, we write the reaction probability for photon exchange in phase- 
space. Comparing the full reaction probability with the reaction rate density for absorbing a photon, we identify 
the effective phase-space photon distribution. This photon distribution is the effective photon number density in 
phase-space and it has the form of a phase-space source folded with a phase-space propagator. We calculate the 
photon number density surrounding a classical point charge and explain how the photon's phase-space propagator 
and phase-space source work. Finally, we comment on the implications of this section for the parton model. We will 
find that we understand how partons propagate and have an idea how to make the gluon distribution gauge invariant, 
but since our photon source is point-like we do not learn anything about the parton sources. 

To find the Parton Distribution Functions, one can solve the parton evolution equations or equivalently one can 
sum up a class of parton ladder diagrams. The simplest parton ladder has one rung corresponding to a single partonic 
splitting. In Section III, we study the QED analog of this process: a virtual photon splitting into a virtual electron 
and on-shell positron. We start our analysis by generalizing the phase-space Weizsacker- Williams Approximation to 
include electrons and writing down the effective electron distribution. This effective electron distribution takes the 
"source-propagator" form. While this "partonic" splitting leads to a complicated form of the electron source, the 
shape of the source is mostly determined by the underlying photon (the "parent parton" ) distribution. We calculate 
the electron distribution explicitly for a classical point charge and discuss how the electron propagates from the source 
to the observation point. We have another reason for studying electrons: we can discuss a simple case where the 
"source-propagator" picture breaks down. This is the case of two virtual photons colliding to produce an electron 
positron pair. The "source-propagator" picture breaks down because the photon fields interfere on the length scale 
of the electron-positron creation region. Nevertheless, discussing the process in phase-space gives us insight into the 
reaction dynamics. Finally, we comment on the implications of the results from this section for QCD parton densities. 

In Feynman perturbation theory, both the photon and electron phase-space densities have a "source-propagator" 
form. This form does not usually arise when one uses time-ordered nonequilibrium methods because one usually 
derives transport theory only after making the gradient approximation. The gradient approximation amounts to 
ignoring small-scale structure of the particle phase-space densities, resulting in much simpler collision integrals [^,^. 
In Section LV, we follow essentially the standard semiclassical transport equation derivation, but never make the 
gradient approximation. Thus, we arrive at the Generalized Fluctuation-Dissipation Theorem which codifies the 
"source-propagator" picture of the particle densities. Crucial inputs to the theorem are the phase-space sources; 
we will discuss how to calculate them. With the sources and the Generalized Fluctuation-Dissipation Theorem, we 
derive a set of phase-space evolution equations. These evolution equations describe the evolution of the system in 
phase-space from the distant past to the present, including all "partonic" splittings, recombinations and scatterings. 



■^Strictly speaking, neither the phase-space sources nor propagators are probabihties as they can be negative. As with any 
other Wigner transformed quantities, they must be smoothed over small phase-space volumes to render them positive definite. 
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Furthermore, we can expand these evolution equations to get the lowest order contributions to the particle densities 
or we can differentiate the evolution equations to get transport equations. 

As a practical application of this study, in Section ^ we examine the coordinate space structure of the parton cloud 
of a nucleon. In principle, one should Wigner transform the quark or gluon wavefunctions of a nucleon. Since we do 
not know the quark or gluon wavefunctions of a nucleon, such a specification is not possible and we must result to 
model-building. One might envision constructing a model phase-space parton density of a nucleon by multiplying 
the momentum space density (the Parton Distribution Function) and the coordinate space density of the partons pl| . 
This approximation neglects correlations between the momentum and position in the parton density which are present 
in the phase-space density |^,^. One might insert these correlations using uncertainty principle based arguments 
|pO| , pT| . This has intuitive appeal, but such a prescription is ad-hoc at best. We can approach this problem in a more 
systematic manner using some physical insight from the momentum-space renormalization-group improved parton 
model. In this model, the parton densities are calculated by evolving the parton densities in virtuality (Q^) and in 
longitudinal momentum fraction (x) . This evolution is equivalent to evaluating a certain class of ladder diagrams and 
these diagrams can be re-cast in the form of the phase-space Generalized Fluctuation-Dissipation Theorem. Thus, we 
can discuss the parton phase-space densities of an hadron in the large-Q^ limit or in the small-x limit. We find that 
neither large-Q^ partons nor small-cc partons extend beyond the nucleon bag in the transverse direction. However, we 
find that the large-Q^ partons extend out an additional^ Hc/xPl from the bag surface in the longitudinal direction. 



This is in line with what others have estimated [gO 21 1. Furthermore, we estimate that the small-x partons extend at 
least an additional 'hc\J —(p- from the bag so the small- a; parton cloud is substantially larger than the large-Q^ cloud. 

Throughout this paper we use natural units (ft = c = 1) when convenient, but we insert factors of he whenever 
directly comparing a length to an inverse momentum. The signature of the metric tensor is (+,—,—,—). 

II. THE PHASE-SPACE PHOTON DENSITY 

One calculates a parton-model cross-section by folding a Parton Distribution Function (PDF) with the cross section 
for the partonic subprocess. One follows a similar procedure for calculating the cross-section in the Weizsacker- 
Williams Approximation |H-|lql: one folds the effective photon distribution with the cross section for absorbing 
the photon (shown in Fig^~0(b)) to obtain the full cross section (shown in Fig. 0(a)). We recast the Weizsacker- 
Williams Approximation in phase-space and, in the process, define the phase-space effective photon distribution. 
This phase-space photon density has the form of a phase-space source convoluted with a phase-space propagator. 

Let us outline this section. First we will compute the photon/current B reaction rate in phase-space. This is 
simply the probability for the probe particle, B, to interact with a free photon. This calculation is simple so we 
use it to illustrate how we rewrite everything in phase-space. Second, we will calculate the reaction probability for 
one-photon exchange. In the Weizsacker-Williams Approximation, the reaction rate is supposed to be the effective 
photon distribution folded with the photon/current B reaction rate, so we can identify the the phase-space effective 
photon distribution. The effective photon distribution is a gauge-independent effective number density of photons. 
Next, we calculate the phase-space photon density surrounding a classical point charge. This calculation will highlight 
how the photon source and the propagator function in phase-space. Finally, we conclude this section with a brief 
discussion of the implications for a phase-space version of the parton model. 

In all of our calculations, we find the reaction probability. Mapping our results to cross sections is trivial and is 
outlined in Appendix Since we are finding reaction probabilities, we work in Feynman perturbation theory. 

A. Photon/Current B Reaction Rate 

We start this subsection by finding the photon/current B reaction rate, W-yB-^B' (a^,?)- This reaction rate plays 
the role of the "partonic" subprocess cross section in a parton model calculation. Our derivation demonstrates how 
to rewrite the reaction probability completely in terms of phase-space quantities. The high point in this calculation 
occurs in equation ( ^.l[ ) when we identify the Wigner transforms of B's current and of the photon field. This type of 
identification lets us rewrite the reaction probabilities in phase-space. 

To find W-yB^B' (2:, 9) we write the S-matrix for the process in wL 



^The nucleon has 4-momentum — {Pq, Pl, Or)- 
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S.B^B' = J d^x (0| A^(.t) |g,A) {B'\3,{x) \B) 

= / d'x^e-^^-^{Q\A^{x)\q,\){B'\j,{k)\B). 
J (,2'"'] 

Here (0| \q, A) — y^ ^|^^^ £^'(A) e*'''^ is the free photon wave function (with = 0) and is the current operator 

for the probe particle B. We leave both the initial and final states of B unspecified so the final state may be a single 
particle or several particles (as in fig. ||(b)). 

We now square the S-matrix and average over photon polarizations: 



\SjB-.B'\ ^ I d Xd X ——^——^e 



I ^ (0| A^(x) \q, A) {q, A| A-(x') |0) {B'\j,ik) \B) {B\ jl{k') \B') 



^2 

A=± 

On writing the coordinates and momenta in terms of the relative and average quantities (i.e. k ^ k — k' and 
K = i(fc + fc')), and taking advantage of the momentum conserving delta functions in the current matrix elements, 

I'S'-yB^B'l^ becomes 

I^.B^ef = / d^Xd^i ^ ^ (01 A^(X + 5/2) \q,X) 

J [ZTr) {Ztt) z [ZA) 

X {q, A| A*-iX ~ i/2) |0) {B'\ ]^{K + fc/2) \B) {B\ jUK - k/2) \B') . 

There are two Wigner transforms in this equation: the Wigner transform of the photon field (the x integral) and the 
Wigner transform of B's current (the k integral). 

Now we rewrite the S-matrix in terms of the phase-space quantities and define the reaction rate density: 



'7B- 



/ 7?T4 ^ E 'f^W <(A) i27rr S\q ~ k) jr(x, k) 

J (27r) Vm\ 



(2.tV\q,\f^^ .w. . , ... , ^^^^^ 

d'^x W^B^B' {x, q) . 
We also have defined the Wigner transform of the current: 

jr(^-,9)^ / ^^e-^'^-^ {B'\f{q + q/2)\B) {B\j^%q- ql2)\B') . (2.3) 
J (27r) 

Since B's Wigner current isproportional to the reaction rate, it is natural to give them the same physical interpretation: 
as a "probability" density,Q| for absorbing a free photon with momentum q at space-time point x. Now, it may not be 
clear where the spatial structure of the reaction rate comes from, especially since the incident photon is completely 
delocalized in space (it is in a momentum eigenstate). To give the reaction rate spatial structure, we must localize 
either the initial or final states of B with a wavepacket. 

B. Photon Exchange 

In this section, we write the reaction rate for one-photon exchange (see Fig. |^(a)) in phase-space. We do it two 
different ways: in terms of the Wigner transforms of the currents A and B and the photon propagator and in terms 
of the Wigner transform of the photon vector potential. The first form of the reaction rate has a clear physical 
interpretation in terms of photon emission, propagation, and absorption. However, it is the second form which can 
be brought into the form of a "partonic" cross-section. 



^Because the Wigner current is the Wigner transform of a quantum object, it may not be positive definite |^,^ so it can 
not be strictly interpreted as a probability. 
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The S-matrix for Fig. |i|(a) is 

Sab-,A'B' = J d\d'y (A'l j^^(x) \A) D;^{x, y) {B'\j^''{y) \B) . (2.4) 

Taking the absolute square of this S-matrix and rewriting it in terms of Wigner transformed currents and propagators, 
we find 



|5'ab^A'B'|^ = \d'^yd'^x^—^j'^{y,q)D'l^^,^,{y-x,q)j!^''{x,q). (2.5) 



Here, the Wigner transform of the photon propagator is 



Dl.^'.'i^^l) = d^xe^^-^Dl,{x + x 12)01* Ax -x/2) 



= (47r)2g^^,g^/^/G"=(a;,g) 

and G'^{x,q) is the Wigner transform of the scalar propagator. We derive G'^{x,q) in Appendix |^. We discuss the 
Wigner transforms of the propagator and current when we study the photon and electron distributions of a point 
charge. 

Equation ( |2.5| ) has an obvious physical meaning: 1) current A makes a photon with momentum q at space-time 
point y, 2) the photon propagates from y to x with momentum q and 3) current B absorbs the photon at space-time 
point X. The spatial structure of the integrand of (^|^) comes from localizing either A or B. 

Now we take a detour and calculate the Wigner transform of the vector potential of the current A. In terms of the 
current density and propagator, the vector potential is^ [ p^ : 



A'^{x)= I d^yDl^[x-y)JAy). (2.6) 

The Wigner transform of this is: 

A,Ax.,q) = i d^ie'^-'>A^{x + x/2)Al{x-S:/2) 



d^y J A " (2/' Dl,^,^,^{x - y, q) 



(2.7) 

, , , '7' — 1/ n\ 



The Wigner transform of the vector potential has a "source-propagator" form. Current A (the photon source) creates 
the photon with momentum q at space-time point y and the propagator takes the photon from y to x. Let us put 



this in equation (2.5) 



\Sab^A'B'?^ j d'x-0^^A^,{x,q)J'^\x,q). (2.8) 

Stated this way, the spatial structure of the integrand of this equation comes from either localizing B or from the 
spatial structure in the Wigner transform of the photon vector potential. 

Equation ( p.Sf ) is close to the form of a cross section in the parton model because current B is proportional to the 
photon/current B reaction rate (as discussed in the previous subsection) and the vector potential is proportional to 
the phase-space effective photon density (as discussed in the next subsection). 

C. The Weizsacker— Williams Approximation 

The effective photon distribution we derive here is the phase-space analog of Weizsacker- Williams' effective photon 
distribution. Thus, it has the interpretation as the number of photons with virtuality q^ in a unit cell of phase-space. 
One could use it to calculate the reaction probability in Fig. |^(a) by folding it with the photon/target reaction rate 
W-yB^B'{x,q). In parton model terms, we take the Parton Distribution Function (the effective photon distribution) 



^Jackson actually uses the retarded propagator to define the vector potential because he discusses classical fields. 
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and fold it with the partonic subprocess (the photon/target reaction rate) to get the reaction probabiUty for the entire 
process. 

We win derive the Weizsacker-Wihiams Approximation in several stages. First we decompose B's Wigner current 
into photon polarization vectors, allowing us to rewrite W-yB^B' {x, q) in terms of J^'^{x, q). Knowing this, we identify 
the effective photon distribution. In the final subsection we will discuss the gauge independence of our effective photon 
distribution. 



1. Current Decomposition 



If the photon probing Jq'^{x, q) is sufficiently delocalized in space (i.e. daAf^,j{x, q) <C qaAf^,j{x, q)), the momentum- 
space cutting rules tell us that we can expand J^'^ {x, q) in terms of the photon polarization vectors [p4[ : 

jr(^,g) = 5]e^(A)e-(A) J,™„.(x,g) 

+ '^e^t^)e*''{0)Jscaiar{x,q) (2.9) 

+ —^Jlong[X,q). 

Here, £^(0) is the scalar (i.e. time-like) polarization vector: €^,{0) — psp, — qfj.q -ps/q^, where ps is the momentum of 
B. The transverse polarization vectors, e^(±), span the hyperplane perpendicular to e^t(O) and g^. Now, if Af^i,(x,q) 
is not delocalized, then Eq. (|]^) should be modified to include gradient^ in x. However, if we were to include those 

gradients here, we could not map J^'^{x,q) to W-ys^B'- 

Since e''(A) e* (A') = Sxy, it is simple to find the separate currents in ( |2.9| ) in terms of J^'^{x, q): 

Jscalar{x,q) = 6^ (0) 6* (0) J^" (x, g) 

and 

X=± 

The longitudinal piece, Jiong{x, q), vanishes due to current conservation. 



The Effective Photon Distribution 



If we insert (2.9) into equation (2.5), the reaction probability is a sum of two terms: 



I^a: 



B-*A' 



+ 



d^q 



A^^x, q) e''(A) e*''(A) Jtrans{x, q) 



(2.10) 



A^,{x,q) e^'{Q)e*''{0)Jscalar{x,q). 



The two terms in ( 2.10| ) describe transverse and scalar photon exchange between currents A and B, respectively. 

Noting that if Jtrans(x, q) has a weak q^ dependence,^ then Jtrans{x, q) oc W^b^b' (a;, q). In other word s, Jtr ansix, q) 
is proportional to the reaction rate for the "partonic" subprocess. Therefore, the transverse term of ( 2.10| ) can be 
written as 



I'S'ab^ 



A'B' 



47r 



4 Vd^qdq^ dny{x,q) 

^"(2^^ d3,rf3grf,2 >^7B^B>(^,g), 



(2.11) 



^These gradients come from Wigner transforming terms proportional to the relative photon mome ntum 
'^The reader should note here that the reaction rates are for photons with any g^, while in Section 
for on-shell photons only. 
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the reaction rate was 
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provided we identify the transverse effective photon distribution as 



dn~^(x, q) 
cPx (Pq dq^ 



(2.12) 



We can make a similar identification with the scalar term. This effective photon distribution is the spin summed 
photon Wigner function. In other words it is the phase-space number density of photons at time per unit q^ . This 
effective photon distribution is the QED analog of the phase-space Parton Distribution Function. This generalizes 
the Weizsacker-Williams method to phase-space. 

While classical derivations of the Weizsacker-Williams method begin with finding the photon power spectra from 
the Poynting flux ||l^,|ljl , a quantum mecha nical derivation follows along the lines of what we do here |^,|6| . Were 
we to perform the spatial integrals in (2.1C), we would find that the exponentials in the Wigner transforms conspire 
to make several delta functions. The resulting delta function integrations are trivial and we would quickly recover the 
momentum-space result. 

Now, multiplying the photon phase-space density by the projection tenso r £^'(A) e*^{X) in ( 2.12 ) does not 

render the photon distribution gauge invariant, unlike in momentum space [ p4[ , because the photon distribution is 
not completely delocalized. When A^u{x, q) undergoes a gauge transformation, terms proportional to are removed 
but terms proportional to d/dx'^ are n ot. A gauge invariant virtual photon distribution is introduced below. This 



gauge invariant distribution reduces to (2.12) when A^jy(a;, g) is sufficiently delocalized 



3. Gauge Issues 



Parton densities are supposed to be gauge invariant but our effective photon distribution is gauge dependent. In this 
section, we discuss how A^^(x, q) transforms under a change of gauge, determine the gauge invariant part of A^i^(x, g), 
and state how the gauge invariant part of A^i,(x,q) is related to the effective photon phase-space distribution. 

In the energy-momentum representation, gauge transforming the photon field adds an arbitrary function in the 
direction of the photon momentum to the photon vector potential: Af^{q) — > Af^{q) + q^if{q)- Because components 
of Af^{q) in the direction of q^ are gauge dependent, we can write Afj^{q) as a sum of the gauge independent and 
dependent parts: 

A,{q) ^ Aliq) + A^^{q) 

where ^ii(q) = ^^A'^(g) is the gauge dependent part of A^{q) and A^{q) = A^{q) — A}li,{q) is the gauge independent 
part. Wigner transforming the photon field gives us a term that is gauge independent and terms which are gauge 
dependent: 



At,,y{x,q) = 



d^q 



Ahq + q/2) +AUq + q/2) 4(g - q/2) + A^ {q - q/2) 



(2.13) 



= Aj^^ix, q) + A^J (x, q) + Ali{x, q) + ^^(a;, q) 



The only gauge independent piece of Ay^i,{x,q) is A-^^{x,q). We do the integrals in (|2.13|) and identify the tensor 
that projects off the gauge dependent part of Ac,p{x^ q): 



= V^,,„pA''P{x,q) 



(2.14) 



where 



{q±id/2)^{q±idl2), 
{q±id/2Y 



(2.15) 



This projector must be understood as a series in q^ and 9^, so can only really be used when So- A^^(a;, q) < qfjAi^i_^(x, q). 
Now, the statement of current conservation for a general J^i/(j:,(?) is 



(g ± ld/2YJ^,{x, q)^{q± id/2YJ^,{x, q) = 0. 



(2.16) 



So, as expected, current conservation ensures that only the gauge independent part of yl^^(.'E, q) appears in the reaction 
probability. 
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With A-~j^(x, q) in hand, we can postulate the gauge invariant photon distribution: 



(2J7) 



This reduces to ( 2.12 ) if the photon field varies slowly in spa ce (i.e. we neglect the g radie nts dcrAf^^{x, q) <C 
QcrAfj_^{x,q)), as we now show. Neglecting the derivatives in ( 2.15D , the projection tensor in ( ^.14 ) reduces to 



9up 



(2.18) 



g^^. Putting 



Since the polarization vectors form a complete basis in Minkowski space, i.e. X)a=± o ^mI'^) ^*u{^)' 
( ^.18 ) in equation ( 2.17 ), we arrive back at ( ^.12 ). 

The tactic of projecting out the gauge dependent parts of the photon distribution works mainly because of the 
simple form of the U{1) gauge transformation. Nevertheless, a variant of this technique probably could be applied to 
gluons. 



D. Photon Phase— Space Density of Classical Point Charge 



We now calculate the density for the simple case of a classical point charge radiating photons. If we localize the 
source's wavepacket and view it on a length scale larger than its localization scale, we can treat its density as a 
delta function. Thus, the shape of the photon distribution is determined the photon propagation and we can use 
this calculation to illustrate how partons propagate in phase-space. We explain that, despite our use of Feynman 
perturbation theory, the photon propagates via the Wigner transform of the retarded (time-ordered) propagator. As 
such, the photon propagates a distance of roughly R\\ = l/\qL\ in the direction parallel to the photon 3-momentum 
and Rj_ = 1/ ^^/\q'^\ in the direction perpendicular to the photon 3-momentum. We demonstrate this behavior by 
plotting the coordinate-space distribution of photons with q^ <C gg (making the photons coUinear with the source) 
and with q^ ~ 9o- Appendix ^ we examine the additional case of a static point charge (i.e. v = 0). This case is 
not relevant for partons as a parton source must be taken in the limit \v\ — *■ c. 



1. Classical Current 



For our source current, we assume the source particle's wavepacket is localized on the length scales that our photons 
can resolve so we can replace its current with the current of a point particle (we discuss when this replacement is 
valid in Appendix |d|) . The source particle follows a classical trajectory a;^ = xijVf^ with four-velocity — (1, vl, Ot), 
vl ~ c and 7 ~ \/{\ — v\) ^ 1. Ignoring the recoil caused by photon emission, the current of the point charge is p^ ] 

ji_,{x) = evi_, S^{x- xqv) . 

The Wigner transform of this is the classical Wigner current: 



^c^iLicaiC^, 9) - y d^S: e^^-^j,{x + x/2)jUx - i/2) ^2.19) 
= 27rae„i Vf^Vi, S{q ■ v)d^{x- Xqv) . 

Here = ttem is the QED coupling constant. 

The current has several easy to interpret features. The first delta function sets q ■ v — 0. This ensures that the 
emitted photons are space-like and that current is conserved. It also insures that, when q^ 0, the photons become 
collinear with the emitting particle {qo = VLqL ~ Ql making « q|i « 0). This delta function arises because we 
neglect the recoil of the point charge as it emits a photon. The second delta function insures that the source is 
point-like and follows its classical trajectory. 

This source has one other feature of note: it allows for emission of both positive and negative energy photons. In 
the following work, we consider creating only positive energy photons so we insert a factor of 26{qo) in ( 2.19|) . This 
amounts to constraining the source's initial energy to be greater its final energy. 
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2. Retarded Propagator 



For our phase-space propagator, we take the Wigner transform of the retarded propagator instead of the Wigner 
transform of the Feynman propagator. This replacement is legal since the reaction probability may be expressed in 
terms of either propagator, provided one uses the same asymptotic states in both cases. Both propagators will give 
different particle densities at intermediate stages but the final state particle densities converge as time goes to ±00. 
We use the retarded propagator here because it leads to a more transparent interpretation fo the particle densities. 
We discuss the Feynman propagator in our discussion of the electron distribution of the point charge. 

The Wigner transform of the retarded propagator gives the weight for a particle with four-momentum to 
propagate across the space-time separation Ax^ = a;^ — y^. In the Lorentz gauge, the retarded photon propagator, 
D^^^/j^, (Ax, q), is proprtional to the retarded scalar propagator G'+(Aa;, q): 



D+^,iAx,q) = g^,g,,,,G+{Ax,q) (2.20) 



The retarded scalar propagator is 



.+ rA^ l^.A^ ^^/A.2^^M2^sin(2VA2) 



G+{Ax,q) -0{Axo)9{Ax^)9{X^)- 



TT 

and it is derived in Appendix In this expression, the Lorentz invariant is A^ = (Ax ■ q)'^ — q^Ax^. 

Let us now estimate how far the retarded propagator can send a particle with the momentum = (qoj'ZLiO)- 
First, the retarded propagator has two theta functions, one that enforces causality and one that forces propagation 
inside the light cone. The rest of the interesting features of the propagator are tied up in the dependence on A^. Since 
G~^{Ax,q) oc 0(A^) sin ('\/A^)/-\/A^, the particle can not propagate farther than the inequalities < < 1 allow. 
To see what these constraints mean, we investigate the q"^ > 0, q"^ < 0, and q"^ = cases separately. 

To study the > case, we position ourselves in the frame where q'^ = {qQ,0). In this frame, the A^ constraint 
translates into a restriction on the spatial distance a particle can propagate: 

< q'o'^Ax'^ < 1. 
Combined with the light-cone constraint, Aaf is constrained to 

I l/kol foi' large Ax'„ > l/\q'„\. ^ ^ 

To find a cutoff for Ax'q, we realize that, for a given g^^, the propagator gives the "probability" distribution for 
propagating across the space-time displacement Ax'^. Thus, we can integrate G~^{Ax' ,q') over all space and over 
time up to some cutoff time r, giving us the total "probability" for propagating to time r. We find that the propagation 
probability becomes unimportant for r > l/|<7ol and this sets a cutoff in Axq: 

Ax', < l/\q'o\. (2.21b) 

Together, these three constraints define the space-time region where the particle can propagate. When we move back 



to the frame with = (go, gi, Ot), the region contracts in the temporal and longitudinal directions. From Eq. (2.21), 
the limits of the propagation region are 

lAfrl < i?± - -4=^ (2.22a) 



|Axl| <i?|| = ^ (2.22b) 

ml 

lAxol <Ro = T^ (2.22c) 

m 

We study the g^ < case in a similar manner. In the frame with g^^ = (0, g^,OT), the A^ constraint implies 

< g'i(Ax'o - Ax%) < 1. 
Combining this with the light-cone constraint immediately gives us a limit on Ax'j^: 
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lAx^l < (2.23a) 

As with the > case, we can integrate the propagator to find the total "probabihty" for propagating to the time 
r. In this case, the propagation probabihty is important only for r < 1/|<Z^|, giving us a limit in Axq of: 

lAx^l < (2.23b) 

The limit on |Aa;^| then follows directly from the light-cone constraint: 

lAf^l < (2.23c) 

Again, these constraints define a the space-time region where the particle can propagate. Boosting back to the frame 
with = (go, ql, Ot), again the longitudinal and temporal spread gets Lorentz contracted: 

\Axt\ < i?± = (2.24a) 
\Axl\<R\\ = j^ (2.24b) 



190 1 

|Aa;o|<i?o = ^. (2.24c) 



Now we study the = case. With = 0, becomes 

= \Ax-q\ = \qo\\Ax ■ q - Axa\ < 1. (2.25) 

On other words, high energy particles tend to follow their classical path while low energy particles can deviate from 
their classical path. Expression (2.25) then gives a measure of the deviation from the classical path. 



3. Phase-Space Effective Photon Distribution 

Now we put these elements together into the photon density. We concentrate our efforts on A^^{x, q) because all of 
the spatial dependence of the photon distribution is tied up in the Wigner transform of th e vector potential. Inserting 
the classical current and the retarded photon propagator in the Lorentz gauge into ( ^.7[ ), we find 

A^,^{x,q) = ^TTaemV^iVvd{qo)5{q- v) j d'^y C^ix - y,q) 5^{y- yov) . (2.26) 

The delta function integrals in equation ( p.26| ) are trivial, however the remaining proper time integral can not be done 
analytically. We find 



A^^(x,q) = i&^r^e,nieiqo)S{q ■ ^) ^^^^^(^^^^ . s^-gV ((^ ' v? ~ x^v^) ~ (^ ' , (2-27) 



(87r)^Q; 

em 

j9{qQ)S{q ■ v) 

where the dimensionless function A{a, b) is given by 



There are two interesting cases that are easy to explore: that of photons nearly coUinear to the source particle (i.e. 
Qa ~ Ql ^ \qt\), and that of photons with a large transverse momentum (i.e. \qx\ ~ qo, qh)- Since there is a g^ 
singularity in the photon density and nearly on-shell photons (i.e. q^ 0) are coUinear, there will be many more 
coUinear photons than any other kind. 

Two plots, representative of coUinear photons and high q^ photons, are shown in Fig. The left is a plot of the 
dimensionless function A{a,b) for coUinear photons with = {me,me/vL,OT)- On the right is a plot of photons 
with transverse momentum comparable to their transverse momentum and energy, q^ — (me, mg/wL, 0.56MeV/c, 0). 
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The characteristic energy scale of QED is rUe, so we choose this scale for the momenta to plot. In both plots, we 
chose vl = 0.9c to illustrate the Lorentz contraction of the distribution. The oscillations exhibited by both photon 
distributions are expected for a Wigner transformed density |2^,^^. To obtain an equivalent classical distribution, 
one should smear this distribution over a unit volume of phase-space. 

Both cuts through the photon distribution show Lorentz contraction. For the collinear photons, this contraction 
occurs in the longitudinal direction. We can account for the contraction with the behavior of the retarded propagator. 
We expect that the width will be ~ = hc/\qQ\ parallel to q and ^ R±^ = hcj \f\cp\ perpendicular to q. For the 
collinear photons, q is in the longitudinal direction and go = 7\/|?I so ^ Rl = hc/^yj\q'^\. In other words, the 
collinear photon distribution is a "Lorentz contracted onion" centered on the moving point source. The inner layers 
of this "onion" correspond to higher \q^ \ photons. However, we must emphasize that the contraction is not due to the 
movement of the source, but rather due to kinematics of the photon's creation and the propagation of the photon. 
To illustrate this point, one only needs to look at the high transverse momentum photons: their distribution is tilted. 
In the case plotted on the right in Fig. |^, the photon momentum points 45° to the longitudinal direction, coinciding 
with the tilt of the distribution . Furthermore, the width of the distribution is ~ R\\ = hc/\qo\ along this tilted axis 
and ^ = he/ yj\q'^\ perpendicular to this tilted axis. 



4. Comment on the Gauge Dependence of the Effective Photon Distribution of a Point Charge 



Now, Afj_^{x, q) is a gauge dependent object and the physically interesting object, the effective photon distribution, 
is gauge invariant. One might ask whether the interesting features of Anjjx,q) disappear under a gaug e transform. 
To see whether this happens, one must inse rt A^ ^lx, q) into equation ( 2.17 ); the only things in ( 2.17 ) that could 
significantly alter shape of the distribution (2.27) are the gradients. Now because the photon source is extremely 
localized (it is a delta function), the shape of the photon distribution comes solely from the propagator. Since the 
propagator varies significantly on length scale comparable to 1/qn , derivatives of A^i/(x,(?) are always comparable in 
size to q^ and any expansion of the gauge projector in Eq. (2.14) will not converge. So, we must conclude that our 
photon distribution can not be made gauge invariant. Now, had we not used a point source for our photons, the 
integration over the source could smooth the photon distribution so that it varies slower in space. In that case, our 
distribution could be rendered gauge invariant. 



E. What the Photons Tell Us about the Partons 



Parton model cross sections can be written in phase-space as a folding of the phase-space parton distribution 
function with the reaction rate for the partonic sub-process. The phase-space parton distribution functions are the 
spatial number density of partons with a certain momentum. The parton distribution functions have a "source- 
propagator" form and can be defined in a gauge invariant manner. If the phase-space parton source produces only 
positive energy partons or if we use time-ordered field theory then the partons propagate from their source using the 
Wigner transform of the retarded propagator. This retarded propagator propagates off-shell partons up to roughly 
~ = /ic/min(|go|, l^l) parallel to the parton three-momentum and ~ R± ~ hcj ^J\cp\ perpendicular to the parton 
3-momentum. Both of these estimates are valid only in frames with go,(f 7^ 0. When either go = or g = 0, 
propagation is cut off at '-^ = ^cj \f\(p\. On-shell (i.e. g^ = 0) partons tend to follow their classical trajectory, 
with deviations from that trajectory of order ~ l/|go|- 

Despite what we have learned, we know next to nothing about parton sources in phase-space. We use a point source 
while a nucleon has spatial structure on the length scales of interest. Furthermore, partons radiate other partons and 
this alters the source. We gain more insight into the phase-space sources in the next few sections. 



III. PHASE-SPACE ELECTRON DENSITY 



In the parton model, the Parton Distribution Functions can be found by summing a class of ladder diagrams and 
the simplest of these has only one rung, corresponding to a single partonic splitting. One can see the QED analog of 
the first rung of such a ladder in Fig. |^(a). Probing the electron distribution occurs in three steps: 1) a virtual photon 
splits into and electron-positron pair with the positron on-shell, 2) the virtual electron propagates from the splitting 
point toward the probe particle and 3) the electron interacts with the probe. We can give this rate a parton-model 
like form by associating steps 1 and 2 with the effective electron distribution and step 3 with the electron/probe 
reaction rate (see Fig. ^(b)). 
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Let us outline this section. In Subsection III A, we demonstrate that the reaction rate for the virtual electron 



exchange process in Fig. y(b) factorizes in phase-space, giving the reaction rate a parton model like form. In the 



process, the electron phase-space density acquires the "source-propagator" form. In Subsection |I1I B| , we calculate 
the electron distribution of a point charge. The electron source shape is determined mostly by the shape of the 
parent photon distribution. The electron can equivalently propagate with the retarded or Feynman phase-space 
propagators. We assume the electrons are massless throughout this section because we only know the form of the 
massive propagators in the high mass limit (me pQ,\p\). This limit is irrelevant for QCD as Aqcd ^ Tnq <^ poJM 



For completeness, we calculate the electron distribution in the high-mass limit in Appendix ^. In Subsection III C , 



we discuss an apparent failure of the "source-propagator" picture: lepton pair production in the strong field produced 
by two point charges. Because the photon fields of the two point charges interfere, it is not possible to clearly isolate 
the source or probe and we can not factorize the square S-matrix into an electron distribution and electron/probe 
interaction. Nevertheless, we can still discuss the process in phase-space, even though we cannot write down the 



electron distribution. Finally, in Subsection IllD we discuss the implications of this section. The discussions of 



splitting, of the massless propagators, and of a failure of factorization are all relevant for partons. 

A. Factorization and the Effective Electron Phase— Space Distribution 

First we show that the process in Fig. |^ can be factorized in phase-space, giving a parton model-like form. The 
S-matrix for the process in Fig. ^(a) is: 



S'^B-toB' = j d xd y A^j,{x)ii)e{x, s)ej^S''{x - y)VBc^B'(2;)- (3.1) 

The spatial structure of the electron source comes from localizing the photon vector potential, An(x). The spatial 
structure of the "partonic" subprocess comes from Vbb^b' (y), the electron/probe interaction in Fig. ||(b). In Eq. (^), 
^g{x, s) = / (0 ^(^' ^y''"^^ is the final positron wavepacket and S%x - y) = J ^ e-'P<^-y \/+J^l^^ is the 
electron Feynman propagator. 

We square iS'^b^cB' and write it in terms of phase— space quantities: 

2 f d'^p d'^q d*k 

I^^B^oB'l = ae„i / d'^x — — 4 — — 4 — 

J {2tt) {2n) {2n) ^ ^32) 

xTr |i(^ - m)rS^{y - x,p)Vbc^b' (y,p)7" 

Here, S'^ijj — x,p) is the Wigner transform of the electron Feynman propagator and it can be written in terms of the 
Wigner transform of the scalar Feynman propagator, G'^{x, q): 

J (zTTj 

= ^ + TOe)a/3(^- « ^ + TUe) a' f3' G" {x, p) . 

Also in Eq. (^^), Vbc^B'(2/,p) is the Wigner transform of the electron/probe interaction and f{x, k) is the phase-space 
density of final state positrons. 

Since the positron rung in Fig. ^ is cut, wc can put the final positron in a momentum eigenstate^ and sum over the 
positron final states and spin. Furthermore, we can separate off the spinor structure of the electron propagator and 
shift the derivatives to act on Vbc^B' ■ In the end we find 



IS'^B^cB'l^ = Otem J d'^xd 



4™ J4 '^^P d'^fc/ 

^ ^4 /r, \4 r, I , I /r, \3 



(2^r (2^) 2|fc/o|(2^)^ 
xA^,ix,q)G^iy-x,p){27:)^S\kf+p~q) (3.3) 

xTr(i(^/ -m)7^(^ + z ^/2 + m,)VBe^B'{y,p){^ ^ i + rueh" 



®This makes the positron momentum weight-function /*(fc) oc S'^{k — kf), with k^ = m^, and the positron phase-space density 
fix, k) = {2V\kso\)-\2T,f5\k - kf). 
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Since Vbo^b' {y,p) is separated from the electron propagator, the reaction probability factorizes. We could explicitly 
calculate the rate for the "partonic" subprocess eB B' , but we are only interested in the shape of the distribution 
as a function of electron momentum. We can guess the form of the electron density just by looking at equation ( [3.3| ), 
without performing the explicit rate density calculation. The electron density is 



d^yd^pdp^ J J 2|fco|(27r)'' (3 4) 

d'^x G'^{y — x,p) Y.{x,p). 



Here, G^{x^p) is the Wigner transform of the scalar propagator. Equation (3.4) has the "source-propagator" form: 
the integral of A^^{x, k+p) over the positron momentum, S, plays the role of the "partonic" source. Because the 
emitted positron is in a momentum eigenstate, the spatial structure of the source comes solely from the parent photon's 
phase-space distribution. 

At this stage, we see several important features of the source. First, we note the d'^k/\ko\ in the positron momentum 
integral. This factor weights positron emission toward small fco. In a typical parton ladder, ordering momenta to 
maximize the contribution from this particular singularity leads to the so-called BFKL evolution equations p5| , p6t . 
Second, we note that the entire spatial dependence of the electron source comes from the parent photon distribution. 
These two points are especially important for Section ^ so they are elaborated on in the next subsection. 



B. The Effective Electron Distribution of a Classical Point Charge 

Our main interest is with how the "parton ladder" (in our case, the source is only one rung of the ladder) shapes 
the electron distribution. First, we discuss the electron's source and how both the parent photon distribution and the 
cut positron rung effect it. Second, we discuss the interplay of the electron creation and propagation. Because the 
electron has positive energy, we can use either the retarded or Feynman phase-space propagator. We choose to use 
the retarded propagator but, for completeness, we describe the Feynman propagator. 



1. The Electron Source 



For our electron source, we choose the photon distribution of Eq. (2.27). This is not a parton-like parent distribution 



as the photon source is poi nt -l ike. Nevertheless, we can use it to illustrate the general features that one expects from 



the electron source in Eq. (3.4). In particular, we discuss the l/|A:o| singularity from the positron rung and we detail 
both the shape of the source and how this shape depends on the photon distribution. 
Up to irrelevant constants, the electron source is 

S(x,p) cx aem / d^k 9{-ko)e{po + ko)d{k^ - mj) 6{q ■ v) (3.5) 
J V-{k+p)^ 



where a — 2\x ■ {k + p)\ and b = 2^J —{k + p)^7^((x • w)^ — x^w^) — (a; • (fc + p))^. The longitudinal and temporal 
positron momentum integrals can be done with the delta functions, leaving the transverse momentum integrals: 



Y.{x,p) (X aemO{p ■ v) 



d^kr 



|fcT|<fcrmax V^Tmax ~ kj 

^ I 0(j)Q + ko+)A{a+,b+) ^ 9{pa + ko-)A{a-,b-) \ 
Here, k^^^^ = 7^(?' ■ v)'^ ^ nil- The two roots of the positron momentum are given by 



(3.6) 



ko± = -7(7P-WTWLV^Tmax-^T) ^3^^ 



kL± = -l{lVLP-vT\/kl^^^-k^). 

Now, in a parton ladder we expect to find a factor of d^k/\ko\ for each cut rung. Here is no exception, one can 
see that d'^k6{—ko)6{k^ — ml) gives us this factor. However, because we neglect the recoil of the source, we have an 
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additional 6{q ■ v) and the factor becomes <Pk/ ^Jk?^^^^ — k^. Because Ikr] < fcrmax = • v <^ \ko±\, I^Lil, this 
singularity forces the positrons to be anti-coUinear with the point charge. Furthermore, because « (because of 
the g2 singularity) the electrons are coUinear with the point charge. 

As in a parton ladder, the shape of S comes froms the parent's distribution. In the case at hand, we can actually 
estimate the (q^) that gives the dominant contribution to S. Because fc^ and kg are fixed by the delta functions 
and kx is bounded by fcTmax, we can estimate the average positron recoil momentum. The average kx is given by 
i ^fcTmax- In general, for vl ~ 1, the average (fcip^) is given by 



(fc^^) ^j\p-v) (^-l±i,-l±i,cos(0T)/7,sin(0T)/7 

For the purposes of illustration, wc choose to emit the positron in the direction kr-xr — cos{9t) — By momentum 
conservation, the dominant photon momentum is (g±^) = p^j. + {k±p). 

On the left in Fig. |[ we plot the electron source for p^ = (2.0, 2.05, Ot) MeV/c electrons from a point charge moving 
to the right with vl = 0.9c. We choose this p^ because it is both collinear with the point charge and because it is 
space-like {p^ < 0). Our source can emit both p^ > and p^ < electrons, however the typical parton in a parton 
ladder is either space-like or on-shell. On the right in Fig. ^, we also plot the photon distribution corresponding to 
the dominant (qp,). Note that both the source and the photon distribution have approximately the same width in 
both the longitudinal and transverse directions. The tilt in the photon distribution gets averaged away in the kx 



integrals in equation (3.6) 



2. Electron Density Using the Retarded Propagator 



Now we put elements of the electron distribution together. In Eq. (3^), we need the Wigner transform of the 
Feynman propagator. However, since the electrons have positive energy w e can r eplace the Feynman propagator 
with the retarded propagator. We discuss the retarded propagator in Section II D 2 and we describe the phase-space 
Feynman propagator in the next subsection. 

We are interested in electrons that have momenta that are both space-like and collinear with the source (for 
comparison with partons), so we plot the coordinate space distribution of electrons with = (2.0, 2.05, Ot) MeV/c 
in Fig. ^. The point source is moving to the right with velocity 0.9c. Both the source and the underlying ph oton 
distribution for these electrons is shown in Fig. H. To perform the four-dimensional spatial integral in Eq. ( |3.6[ ), we 
use a Monte-Carlo integration scheme ||2^ . This integration scheme, being probabilistic by nature, returns both the 
integral at a point and the error on the integral at that point. The nonzero data points never had a relative error 
greater than 20%, but due to this error, the location of the zero contours is uncertain by ^ 30 fm. 

Comparing the electron distribution with the source, we see that the electron distribution is elliptical with longitu- 
dinal and transverse widths compara ble t o what one expects by adding the source width in Fig. to our estimates 
for the propagation length in Eqs. (2.24). Unlike the electron source distribution, the electron distribution is not 
symmetric about xl = 0- This is caused by the positron recoil because, were there no positron recoil, we would have 
a delta function to insure po = PlVl (as we found for the photons). Because of the positron recoil, the delta function 
is widened and the additional spread in energy causes the electron to preferentially propagate forward. 



3. The Feynman Propagator 

Even though we choose situations where we can avoid using the phase-space Feynman propagator, we should 
describe how it works. While the Feynman propagator propagates a particle with a given momentum (say = 
(Po,PL:Ot)) across a space-time displacement Aa;^ — {Axq, Axl, Axt), it does so in a manner very different from 
the retarded propagator. The Feynman propagator is 

G'=(Aa;,p)= ^ [sgn(Ax2) + sgn(p2) + 2 sgn(Aj: • p)] 
The combination of the sign functions in the square brackets can be rewritten in a more transparent form: 
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4 if Ax -p, p^, Az^ > 

-4 if Ax • p, p^, Ax^ < 

2 sgn(Ax • p) if p^, x^ have opposite sign 

Thus, particles with time-like momentum tend to travel forward in time and inside the light -cone and particles with 
space-like momentum tend to travel backwards in time outside the light cone. Also as one might expect, anti-particles 
with time-like momentum tend to travel backwards in time inside the light -cone and anti-particles with space-like 
momentum tend to travel forwards in time outside the light-cone. 

The rest of the interesting features of the Feynman propagator are tied up in the d ependence on the Lorentz 
invariant = (Ax • pY — Ax^p^. As with the retarded propagator in Subsection II D| , we will study the p^ > 0, 

< 0, and = cases separately. 

To study the p^ > case, we boost to the frame where p'^ — (pg, 0). In this frame, = p'q|Ax'P > 0, so only the 



sine term contributes. The sine term is greatest for v A^ < 1 so we have the following limit on the spatial propagation 
distance: 

lAx'l < ^. (3.8a) 

As with the retarded propagator, we can compute the total "probability" to propagate to certain time. This calculation 
gives us the following limit on the temporal propagation distance: 

lAxol < ^. (3.8b) 

Boosting the space-time region defined by these constraints back to the frame with p^ — (po,Pl,Ot), we find the 
following constraints: 

\Axt\<R± = ^= (3.9a) 

|AxL|<i?||=^ (3.9b) 

|Axo| < i?o - A (3.9c) 
iPol 

These limits are exactly the same as the ones we found for the retarded propagator in Subsection |I 

To study the p < case, we boost to the p'^ = (0,p^, Ot) frame. In this frame, A^ = p'^(Ax'o — Ax'^). Inside the 
light-cone, the exponential term disappears and we get a constraint on A^: 

> A^ = p'liAx'l - Ax't) < 1. 

We can integrate to find the total "probability" to propagate to a certain time, giving us a limit on Axq: 

lAx^l < (3.10a) 



Using the A and light -cone constraints, we find similar limits on Ax^ and Ax^: 

1 

1 



Ax^l < ^ (3.10b) 
Ax^l < (3.10c) 



Boosting back to the = {po,pl, Ot) frame, we find 

|p2 



Axt\<R± = ^= (3.11a) 



|Axl| < i?|| = 7^ (3.11b) 



1 

w 



|Axo|<i?o = ^, (3.11c) 
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which is what we found for the retarded propagator. Now, outside of the Hght-cone the situation is more comphcated 
and we must integrate the propagator in the various directions to find hmits. We find: 



(3.12a) 
(3.12b) 
(3.12c) 



When we boost back to the frame with = (po,pj, Ot), we find the result in Eq. ( 3.11c ) 
Finally, we investigate the — case. With = 0, becomes 



< A' 



P\ 



\pq\\Ax ■ p - Axo\ < 1 



(3.13) 



because the exponential term does not contribute on the light cone. On other words the Feynman propagator functions 
exactly like the retarded propagator: high energy particles tend to follow their classical path while low energy particles 
can deviate from their classical path. Expression ( pl3| ) then gives a measure of the deviation from the classical path. 

We find that, despite the different boundary conditions on the two propagators, both the Feynman and retarded 
propagators send particles the same distances. This is probably no surprise since a calculation done using Feynman's 
formulation of perturbation theory must give the same results as the same calculation done using time-ordered 
perturbation theory. 



C. Failure of Factorization and the "Source— Propagator" Picture 



In this section, we investigate electron-positron pair production in the strong field of two point charges. One 
might visualize this interaction as a virtual photon from one point charge probing the virtual electron distribution 
of another point charge. Thus, the electron distribution would appear factorized from the virtual electron-virtual 
photon collision process. However, we will show that this picture is incorrect because the photon fields interfere with 
one another on length scales comparable to the size of pair production region. Of course, this also means that our 
"source-propagator" picture fails here. Nevertheless, we can still formulate the problem in phase space and discuss 
the interplay of the interaction length and particle production length scales. 



1. Interference of Photon Fields 

We can write down the S-matrix corresponding to the process in Figs. ^ and |^ using the same procedures used in 
the previous sections. To lowest order in the coupling strength, we obtain: 

J (27r) (27r) (27r) 

(3.14) 



xA^^(fci, si, fc2, S2,p) {Al{xi)A''2{x2) + A':^{xi)A'i{x2)} . 

Here xi and X2 are the interaction points of the photons and should not be confused with the classical source particles 
1 and 2. We have already separated the 7766 effective vertex 

Ap,.(fci,si,A;2,S2,p) = M(fci, si)7piS"=(p)7t.w(fc2, 52)- 

In A^t^(fci, si, k2^ S2,p), S'^{p) is the momentum-space Feynman electron propagator. The final state electron-positron 
wavepacket is /*(fci, k2) and we will assume the final ee pair to be free and use t he fre e wavepacket from Appendix 
^ The reader should note that we ca n already see the photons interfering in Eq. (3.14). 

As usual, we can rewrite Eq. ( 3.14 ) in terms of Wigner transformed quantities. However, due to the photon fields 
interfering, the structure of the cross terms are complicated. The |5'i2^i'2'ee|^ is: 
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xf{R~r/2,h,R + r/2,k2) 
xAfj,fj,>„„^{ki,k2,p,r){2Tr)'^S'^{qi + 92 - fci - 

X { {2Tr)HHqi ~h+p) Af {R - r/2, gi) A^' (i? + r/2, q^) 

+ (271)4^4(5^ -~k2-p) Af [R + r/2, q,) Af {R - r/2, q^) (3.15) 
+ y d^f exp [ir ■ {-p H ^—^ — -) - ir ■ {qi - q2)] 

X A''/{R-r/i,qi)Af{R + P/i,q2) 

d'^f exp [ir ■ {-p + ^ ^ ^ ) + ir • (gi - 92)] 

X Af' (i? + r74, gi) Af {R ~ r74, 92) } 

This equation could look simpler if, in the interference terms, we Wigner transformed Ai together with A2. However 
then we would have a virtual electron being emitted by some interference field and then reabsorbed by another 
inter feren ce field and the resulting equations would be impossible to interpret using our photon distributions. In 
Eq. (3.15), we neglect k relative to k in the effecti ve ver tex and in the factors of {2koV) because the final state wave 
packets are sharply peaked in momentum. In Eq. (3.15), R is the center of the interaction points xi and X2 and r is 
the space-time separation of these points. The final state Wigner density is 



1 1 f d^ki d^k2 

{2ki,QV) {2k2,oV) J (2^ (2^ 
x/*(fci + fci/2, k2 + fc2/2)/(fci - fci/2, k2 - k2/2) 

and the Wigner transform of the effective vertex is 

e'P"^A^^(fci, A:2,p + p/2)A;,^,(fci, fc2,j5 - p/2). 

[2n) 

We can write the effective vertex in terms of the scalar Feynman propagator, 

A^^'i.i.'(fci, fc2,p, r) = u{ki, si)jf^{^ + ^ ^ + me)j,yv{k2, S2) 

xw(fc2,S2)7i^'(l*- I ^ + me)jf,'u{ki,si)G''{r,p) (3.16) 
= A^^'i.^'(fci,fc2,p,r)G"=(r,p) . 

We simplify the reaction probability by summing over the final state electron and positron spins. We simplify things 
even further by working in the ultrarelativistic limit, namely when « u| « 0. Under these approximations, we find 

le |2 2 /"j4dj4 d^ki d'^k2 d^qi d'^q2 d^p 

="™y''''^(^(^(^(^(^ 

x/(i?-r/2,A:i,i? + r/2,/c2) 

X ^ \t,t,'pv'{ki,k2,p,r)G''{r,p) (27r)''(5'*(qi + q2 - ki - k2) 

spins 

X < rZTT\'0'irH ~ IfA + T)) /\'^'~ i H ~ r rz. fh \ i H + r rz.n-^] (3- 17) 



{ {2TTf5\qi ~ki+p) Af {R - r/2, qi) Af {R + r/2, 92) 
+ {2nfS''iq, ~k2-p) Af {R + r/2, q,) A^ {R - r/2, 92) 



d^f cos [r ■ {~p + -i-y— ^) - f • (91 - 92)] 
Af{R^flA,q^)Af{R + flA,q2)]. 



X 



Given the relatively simple form of this equation, one would think that we could identify the exchanged electron's 
phase-space density. In fact, if we use free particle distributions for the fin al st ate electron and positron and sum 
over final states, we can identify the virtual electron distribution (equation ( |3.4|) ) in the direct terms. However, we 
can not make the same identification in the interference term and factorization is not possible here. We might find 
factorization again if we had several point charges as one can envision a situation with many photon sources screening 
the photons (a plasma for instance). The photon field might then be an incoherent superposition of photon fields. In 
the absence of photon interference, we might be able to define an effective electron distribution. 
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2. The ee Production Region vs. the Interaction Region 



With equation (3.17), we can discuss the various length scales of the problem. First the ee production region is set 
by the shape and size of the photon distributions. Second, the two photon interaction region's size depends on the 
mass and virtuality of the exchanged electron. 



First, take the virtual photon distribution of the classical point charge from Section [IE. Now, the lowest energy 
and momentum that each of the interacting photons can have i^ q = (me, me/i'L, Ot). Because the high energy or 
far off-shell photons are closer to the point charge then their lower energy and nearly on-shell cousins, photons with 
the minimum have the largest distributions. So, the geometrical overlap of the high energy portions of the virtual 
photon distribution sets the size of the ee production region. In Fig. |^ we illustrate this: the two ellipses represent 
the edge of the photon distribution and the shaded region is the region where the ee pairs can be created. 

Now, the size of the two photon interaction itself is determined by how far the exchanged electron can travel between 
the vertices in Fig. |7[ For this, we look at the phase-space electron propagator. Assuming massive electrons,Q we use 
Remler's causal propagator. Here the phase-space "probability" for propagating between two space-time points drops 
like e"^™""^ for space-like electrons and like sin2r(-\/p^± TOe)/(-y/p^± me) for time-like electrons. The proper time 
along the electron 4-momentum is r. In the direction transverse to the electron four-momentum, the "probability" 
is zero. Thus, the interaction region has a characteristic length scale of « l/mg. This is comparable to the width 
of the photon distributions, so there is no scale separation. Typically one requires the interaction length scale to be 
much smaller than the characteristic length scale of the particle density in order to justify the gradient expansions 
and allow for a transport description. Because our approach does not rely on the gradient expansions, a transport 
description may still be possible. 

D. What the Electrons Tell Us About the Partons 

In this section, we learned several things about the massless parton phase-space densities. First, owing to the fact 
that the simplest parton ladder contains one rung representing a single partonic splitting, we learned how both the 
parent parton and cut rung affect the parton distribution. The shape of the parent parton distribution determines 
the spatial structure of the parton source. The rung of the parton ladder segment gets cut, putting that parton on 
shell. The integral over final states of this parton is weighted toward giving it a low fco. Second, partons propagate to 
the same distances with the Feynman propagator that we found for the retarded propagator, despite the difference in 
the boundary conditions of the two propagators. Finally, we learned that the "source-propagator" picture of parton 
densities fails when the source particle and probe particle interact, even through quantum interference. Nevertheless, 
we can still discuss the process in phase-space with the phase-space sources and propagators, even if the densities 
have no clear meaning. 

IV. QED TRANSPORT THEORY 

The "source-propagator" picture of particle densities seems both common and physically intuitive. In this sec- 
tion, we see how this picture arises in time-ordered nonequilibrium theory by deriving the Generalized Fluctuation- 
Dissipation Theorem in phase-space. The derivation of this theorem mirrors the steps often used to derive the 
semiclassical transport equations. Namely, we derive the Kadanoff-Baym equations and formally solve them to get 
the Generalized Fluctuation-Dissipation Theorem. Unlike other derivations of the transport equations, we do not 
perform the gradient expansion. Instead, we Wigner transform the Fluctuation-Dissipation Theorem. If we then 
insert the Wigner transformed self-energies into the Fluctuation-Dissipation Theorem, we get a set of phase-space 
evolution equations for the particle densities. These evolution equations describe the complete evolution of the system 
in phase-space from some time in the distant past to the present, including all parton splittings, recombinations and 
scatterings. In principle, the equations are nonperturbative, but we can expand them perturbatively. We demonstrate 
this by recalculating the photon and electron distributions of Sections || and III. Following this, we derive transport 



equations from the phase-space evolution equations. Finally, we describe how the results of this section can be applied 
to parton transport. 



®The distribution of photons with q = (me, m^/vL, Ot) is shown in Fig. g. 

^"The me = case is uninteresting because the two photon interaction always extends over the entire ee production region. 
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For those familiar with the common step s in d eriving se miclas sical transport equations from the KadanofF-Baym 
equations, we suggest skipping past Section IV B to Section IV C . 



A. Green's Functions 

Our derivations begin with the contour Green's function which we define as 



for scalar particles, 



for photons and 



iG{x,y)^lf^{xW{y)) (4.1a) 



iD^,{x,v) = (fA^{x)A,{y)) - (^(.t)) (4.1b) 
iSapix^y) = (T'4!a{x)'4)p{y)^ (4.1c) 



for fermions. The (...) = Tr(p. . .)/Tr(/3) is a trace over the system's density matrix, specified at time Iq ~oo. 
The field operators are taken in the Heisenberg picture. T denotes ordering along the contour shown in Fig. ^ This 
ordering can be written as 

fA(x)B{y) = 0{xo,yo)A{x)B{y) ± 0{yo, xo)B{y)A{x). (4.2) 

The upper sign refers to bosons and the lower sign to fermions. The contour thcta function is defined as 

„ , , j 1 if xq is later on the contour than yo 
&{xo,yo) = I otherwise 

In addition to the contour Green's functions ( |4.1a| - |4Tc ), we define the > and < Green's functions: 



zG>(x,2/) = 


{kxWiy)) 








(4(x)i.(y)) - 




{My) 




(^^a{x)i!p{y)^ 






iG<{x,y) = 


(^*{y)Hx)) 








(i.(y)4(x)) - 




{My) 


iS^pix.y) = 


- Gl'(}{y)i^a{x)\ 







(4.3a) 
(4.3b) 
(4.3c) 
(4.3d) 
(4.3e) 
(4.3f) 

These Green's functions are hermitian and contain the complete single-particle information of the system. For 
example, setting x = y gives us the single particle density matrix. Furthermore, Wigner transforming in the relative 
coordinate, we find the off-mass shell generalization of the Wigner function for the particles: 

f{x,p)^iG<ix,p)^ J d\x~y) e*(^-«)-W<(x,2/) = j d\x - y) e'^^-y^'P {$* (yMx)) . 

We identify f(x,p) with the number density of particle (or antiparticles) per unit volume in phase-space per unit 
invariant mass at time xq: 

dn{x,p) 

f{x,p) - 



d^x d^p dp^ 

The off-shell Wigner function is related to the conventional Wigner function, fo{x,p), through the invariant mass 
integration: 
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dn{xo,x,p) f 2 

In terms of the ^ Green's functions, the contour Green's function can be written as 

G{x,y) ^ 9ix„,y„)G>ix,y) + 9{yo,x„)G<{x,y) (4.4) 

for both fermions and bosons. Furthermore, because of the equal time commutation relations, G^{x,y) — G^{y,x) 
and G{x,y) — G{y,x), for both fermions and bosons. 

We define several auxiliary Green's functions in terms of the > and < Green's functions: the retarded and advanced 
Green's functions and the Feynman and anti-Feynman propagators. We write only the equations for the scalar 
particles. For the retarded and advanced propagators, we have 

G^{x,y) = ±e{±{xo-yo)){G>{x,y) - G<{x,y)). (4.5) 

For the Feynman and anti-Feynman propagators, we have: 

G%x,y) ^0{xo~yo)G>{x,y) + 9{yo-xo)G<{x,y), (4.6a) 
G^ix, y) = eiyo - xo)G> (x, y) + e{x^ ~ yo)G< [x, y) . (4.6b) 

One can also obtain these Feynman and anti-Feynman propagators by restricting the arguments of the contour 
propagators to be on one side of the contour in Fig. p|. 



B. Conventional Transport Theory 

In this subsection, we follow the standard derivation of the transport equations up to the point where we find the 
Generalized Fluctuation-Dissipation Theorem. The procedure is as follows: 1) find the Dyson-Schwinger Equations 
for the contour Green's functions, 2) apply the free field equations of motion to get the Kadanoff-Baym equations 
and 3) solve the Kadanoff-Baym equations to get the Generalized Fluctuation-Dissipation Theorem. 



1. Dyson-Schwinger Equations 



The Dyson-Schwinger equations encapsulate all of the nonperturbative effects in the field theory that are possible 
with only two-point functions. We can write the Dyson-Schwinger equations for the photons, electrons and scalars 
using the contour Green's functions: 

(4.7a) 
(4.7b) 
(4.7c) 

In these equations, we represe nt the coo rdinates by their index, i.e. xi — > 1. We present the corresponding diagrams 
in Figs. [lo|(a-c). In equations ( I.7a| - 4.7c ), the non-interacting contour Green's functions have a superscript. 

The self-energies describe all of the branchings and recombinations possible for the photons, electrons and scalars. 
The self-energies are: 



(4.8a) 



0(1,1') =G°(1,1')+ / fi2d3G"(l,2)Q(2,3)G(3,l') 
Jc 

D^^il, 1') = D%il, 1') + J^d2d3 Dl^,{l, 2) n^''''(2, 3)0^3, l') 
^„;3(1,1') =Sl^{l,l')+ [ d2d3 5L'(l,2)S"''3'(2,3)V/3(3,l') 



Q(l,l') = z(eZa ) y^d2d3G(l,3)r;;^42,3,l')i?^.(l,3) 

+i(2zae™^V') / c;2d3(i4G(l,2)r^;';,(2,3,4,l')i?^^'(l,3)Z?,,.(l,4) 
Jc 

+Qmf(1)<5*(1,1') 

n^,(l,l') = -i(-ie(7^)„0) / d2d3 3^^,(1, 2)T^'J^U2,3,1')S0' (3(3,1) 
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+i[eZ df,) J d2 d3 G(l, 2) r^00,.(2, 3, l')G(3, 1) (4.8b) 

+i{2iaemZ^g^^,) [ d2d3d4G(l,2)G(3,l)r!;;^^,(2,3,4,l')i^^v'(l,4) 
Jc 

+nMF(i)5p.^^(i,i') 

S,/3(l,l') = *(-*e(7^)«a') / d2d3 5„-^-(l,2)r^l^'^(2,3,l')i?p.(l,3) (4.8c) 

Jc 

+ SMF(l)<5a/3(54(l, 1') 

In Figs. |ri|(a-c) , we show all the diagrams corresponding to the non-mean-field terms in equations ( 4.8a - 4.8c ). We 
define the contour delta function S^{x,y) by 

S^{x ~ y) for .To, yo on the upper branch 
(5*(a;,y) = ^ for xq, yo on different branches 

-6'^{x — y) for xq, j/o on the lower branch 



2. Kadanoff-Baym Equations 



The free-field contour Green's functions satisfy the equations of motion: 

{dl + M^)G\x,y) =S\x,y) 

{i - me)S'"^(a;, y) = S^fiS'^ix, y) 
Combining these with the Dyson-Schwinger equations, we have 

{df + m2)G(1, 1') = (5^(1, 1') + J^d2 Q{1, 2)G(2, 1') 

dlD^,{l, 1') = Ang^j\l, 1') + 47r^ rf2 n/(l, 2)D,,,{2, l') 

(z ^i-me)^o;5(l,l') =S^p5\lA')+ [ d2i:^t3,{l,2)Sp,f3{2,l'). 



There is a conjugate set of equations (4.£), (4.10) with the differential operators acting on 1' 



(4.9a) 
(4.9b) 
(4.9c) 



(4.10a) 
(4.10b) 
(4.10c) 



Restricting ti and ti' to lie on different sides of the time contour in Fig. we arrive at the Kadanoff-Baym 
equations. 



(92 +M2)G^ (1,1') = f d^X2QMF{xi,X2,ti)G^(x2,tul') 

/■tl 

+ / d2 (Q>(1,2)-Q<(1,2))G^(2,1') 

Jtn 



(4.11a) 



+ rd2Q^(l,2) (G>(2,1')-G<(2,1')) 

Jto 



-^dfD%il,l') = j d^X2nMF{xuS2,h)D%{x2,h-l') 



to. 



d2 n>/'(i,2)-n</'(i,2) i?<,(2,i' 



+ / d2n^/'(l,2)(i^>J2,l')-i?<,.(2,l')) 

Jto 



(4.11b) 
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d2 



d2^ 



to 



(S>^,(1,2)-S<^,(1,2))5|^(2,1') 
1^,(1,2) (5>^(2,1')-5<^(2,1') 



(4.11c) 



Here the > and < self-energies have the same meaning relative to the contour self-energy as the > and < Green's 
functions do relative to the contour Green's functions. Again, there is a set of conjugate equations with the differential 
operators acting on 1'. 



3. Generalized Fluctuation-Dissipation Theorem 

Now we define the retarded and advanced self-energies for scalars: 

0± (1, 2) = QMFS{h,t2) ± e {±{h - t2)) (g> (1, 2) - Q<(1, 2)) 

The photon polarization tensor and electron self-energy are defined in a similar manner. 
Using these, we simplify the Kadanoff-Baym equations: 

poo 

{dl + M^)G^ {1,1') ^ / d2Q+(l,2)G^(2,l') 

/"OO 

+ / d2g^(l,2)G-(2,l') 

J to 

-dfD%{l,V)^ d2n+/'(l,2)i?,\(2,l') 

poo 

+ d2n^/'(l,2)I?-^(2,l') 
(z ^^-m,)S%{l,l') = / d2I]+^, (1,2)4^(2,1') 

J to 

+ d2S|^, (1,2)5^,^(2,1') 

Jto 

If we subtract the > equations from the < equations and multiply the resulting equations by ±^(±(^1 
a second set of differential equations: 

/■OO 

{df + Af2)G± (1, 1') =S\1- 1') + / d2 g± (1, 2)G± (2, 1') 

Jto 

-9?z?±,(i,i') =sHi-i') + d2n±/(l,2)i^X(2,l') 
{i h-^e)S^p{l.l') ^5\1-1')+ rd2S±^,(l,2)5±^(2,l') 

Jto 



Solving the initial value problem posed by equations (4.13) using equations (4.14), we find: 



G«^(l,l') = / d2 d3G+(l,2)Q^(2,3)G-(3,l') 

Jto Jto 



+ J d^X2d^X3 G+il,x2,h)G^{x2,to,x3,ta)G-{x3,to,l') 
D%{1X) ^ d2 d3 7^+^,(1, 2) '''^'(2,3)i^;,,(3,l') 

Jto Jto 

+ [ d^X2 d^X3 D+il,X2,to) D^^^'"' ix2,to,X3,to) D-,^{x3,to,l') 



(4.12) 

(4.13a) 
(4.13b) 
(4.13c) 
(4.13d) 
(4.13e) 
(4.13f) 
ti')), we get 

(4.14a) 
(4.14b) 
(4.14c) 

(4.15a) 
(4.15b) 
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S%ihl')= d2 d3 5+„,(l, 2) (2, 3)5^,^(3,1') (4.15c) 

J to J to 



These equations are the Generahzed Fluctuation-Dissipation Theorem. They describe the evolution of a density 
fluctuation (given by the > and < Green's functions) from Iq to ti. 



C. Phase— Space Generalized Fluctuation— Dissipation Theorem 



We now translate the fluctuation dissipation equations (4.15) into phase-space. We only do so for the scalar 
equation because the photon and electron equations are similar. First we extend the integration region to cover all 
time: 

G^{xi,xi>) = d'^X2 d'^xs G^{xi,X2)Q^{x2,X3)G^{x3,xi') 



+ lim d X2d xsSito- X2o)S{to- X30) 

to— ►-oo J 

xG+{xi,X2) G^{X2,X3) G^{x3,Xy) . 



Next, we Wigner transform in the relative variable xi — xi'i 



G^x,p) = I d^x' ^G+{x,p-x',p')Q^x',p') 



lim I d-'x' -^^G+{x,p;x',p')G^{x'J) 



Xf. — > — OD 



(4.16) 



{2^r 

We recognize the Wigner transforms of the self-energy and initial particle density: 

g^(x,p) = J d'^xe'P-^Q^{x + i/2,x-x/2) (4.17) 

and 

S{to~xo)G^(x,p} = 

d^i e'P-'^S{to - {xo + io/2)) 6 {to - (xq - io/2)) G^ {x + S:/2,x- x/2) . ^^'^^^ 



The delta functions render the initial density independent of pq. We have also defined the retarded propagator in 
phase-space: 

G+ix,p;y,q) ^ [ d^x' d'^y' e'<-P-^'-'^-y'^G+{x + x' /2,y + y' /2)G-ix ~ x' /2,y - y' /2) (4.19) 



At this point, one usually applies the gradient approximation to equation (4.16), eliminating the d'^x' integral. We 
do not do this. 

Next, we assume the translational invariance of the advanced and retarded propagators. This is reasonable at lowest 
order in the coupling since the free field advanced and retarded propagators are translationally invariant. Making this 
approximation, the retarded propagator in phase-space becomes 



G+(x,p;y,q) = i2nfs\p-q) j z e'^'^ G+ {x - y + z /2) {G+ {x - y - z /2)y ^^^0) 
^ {2,,Y5\p-q)G+{x-y,p). 

We will use G'^{x — y,p) in all subsequent calculations. In practice, we will only use the lowest order contribution 
t o G~ ^{x — y,p). This means that we dress the ^ propagators but not the ± propagators when we iterate equation 
(4.15). Thus, our particles propagate as though they are in the vacuum. In Appendix ^ we calculate the lowest order 
contribution to G'^^x ~ y,p). 
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Repeating this for the photons and electrons, we arrive at the phase-space Generahzed Fluctuation-Dissipation 
Theorem: 



G^{x,p) = / d^yG+{x-y,p)Q^{y,p) 



lim 



d'vG+ix-v^p) G^{y,p} 



D%{x,p) = / d''yD+^^,^,{x~y,p)n^^'' {y,p) 



lim 



d'yDl^^,^,{x~y,p)Di^^ {y,l 

d^y s^pa'p- - p) s%f3' {y, p) 



lim 

yo^-ca 



(4.21a) 



(4.21b) 



(4.21c) 



These equations describe the evolution of the particle phase-space densities from y^ — > — oo to the time xq, including 
particle creation and absorption through the particle self-energies. They clearly have the "source-propagator" form, 
but also contain information about the initial particle density. The derivation of these equations does not rely on the 
form of the self-energies and the general form is shown diagrammatically in Fig. |l^. Thus, these equation can be 
re-applied to QCD. We exploit this fact when we discuss the shape of a nucleon's parton cloud. 



D. Phase— Space Evolution Equations 



The first step toward getting the phase-space evolution equations from the Generalized Fluctuation-Dissipation 
Theorem is to calculate the self-energies (i.e. the sources). To do this, we insert equations ( 4.12| ) and (4.4) into the 
self-energy equations and keep only the lowest order approximation to the vertex functions. Thus, we assume that 
the interaction time is much smaller than the other time scales in the problem. So, we arrive at the creation and 
absorption rates: 



g^(l, 1') = laemZ^ di^ G^(l, 1') di'. D^^^-'il, 1') + q|p(1)J4(i _ 1') 
n|,(l, I') = laemTr {-/^S^ih 1') 7.5^(1, l')} + *«e,„^' 9iM G^(1, l') dvu G^{1, l') 

S|^(l, 1') = -zae™(7jaa'5|,^,(l, 1') (7.)^,^I?^^''(1, 1') + j:^p{l)6^^S\l - V) 



(4.22a) 
(4.22b) 

(4.22c) 



We have neglected the second scalar term in the polarization tensor and the second photon term in the scalar self- 
energy because they en ter w ith a factor which is higher order than the other terms we kept. 

The self-energies in (4.22) can be Wigner transformed. Taking care to integrate the derivative scalar couplings by 
parts, we arrive at 



Q^{x,p) 



^%{x,p) 



iaemZ'^ [ j-^ j-^ {qi+q2-i d /'2)f,G^{x,qi) 
J (27r) (27r) 

X (qi + 92 - * a /2),D^^^''ix, q2) {27r)^S\p ~ [q, + q^)) + Q%^{x) 
d'^qi d^q2 



(4.23a) 



Tr 



[j,S^{x,qi)j,S^{x,q2)} (27r)V(p- (gi +92)) 



d^qi d^q2 



{27rf {2n)' 



{qi+q2+i d /2)^G^{x,qi) 



X (<Zi + 92 + « a /2).G< {x, 92) {2n)*S\p - (91 + 92)) + Il<p{x)g^, 
d^qi d^q2 ^> 



,^ ,4 ,^ .S<f,{x,q,){^,),,pDi>^%x,q2) 
(27r) (27r) 



(4.23b) 



(4.23c) 



X {2^f6\p - (91 + 92)) + ^^A^)5c.p 
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We can rewrite these equations direct ly in terms of the particle and antiparticle densities]^ to make their structure 
explicit. In accordance with equations ( 4.3a -f), we define the particle densities as follows: 



iG^{x,p) — {x,p) — scalar densities 



(4.24a) 



iD^j^{x,p) = d^jj ix,p) = photon densities 



(4.24b) 



'q/3 



(x,p) — lepton densities. 



(4.24c) 



Here the positive energy part of the > Green's functions correspond to the density for emission of (|po|j??) quanta, 
while the negative energy part corresponds to the density for absorption of (— |po|,p) quanta. Similarly the positive 
energy part of the < Green's functions correspond to absorbing (|po|iP) quanta and the negative energy part to 
emission of (— |po|iP) quanta. Thus, 0{po)g^ (^tP) is the scalar density and d(—pQ)g^ ix,p) is the antiscalar density. 
We can make simil ar ide ntifica tions for the photon]^ and lepton densities. 

Now, combining ( 4.21 ) and ( 4.23| ) and inserting the particle densities, we arrive at 



9^ {x,p) 



ix,p) 



,4^ d'^qi d'^q2 



,^ ,4 .G+{x-y,p){2^Y5\p-{q,+q2)) 
(27r) (27r) 



Xae^^Z^{ql+q2-^ d I^Tg^ iy,qi)iqi+q2-l 8 /2Yd%{y,q2) 



J d!^yG^{x-y,p)iQ^^{y) 

d^yG+{x-y,p)g^ {y,p) 



- lim 



d%ix,p) = / d'^y 



^ ^ Dt.'A^ - y^P) (2-)'^'b - ('^i + 92)) 
[Ztt) [Ztt) 



X <^ Q!p„,,Tr 



7^' (.y-,qi)Y iy>Q2) 



{ 

+ aernZ^qi + q2 + i 8 /2f g^ {y, qi) (gi + 92 + « 9 /2fg^ (y, 92)} 



+ / d''yD;^^,^,ix-y,p)iU<^iy)g^'' 



lim / d'yD+^,ix-y,p)d^^^'-'{y,p) 



d*qi d'^q2 



d'^y 



(2^)' (2^)' 



SZ0o.'fA^-y^P) {2nf5\p-{q^+q2)) 



^ae,n{Y')a'a"S<npn {y,qi) [Y) fi" f3' d% (^,92) 

+ I d^yS+p^,p,{x-y,p) (±iJ:^-p{y)^ 5^,^, 
d^y 5+^a'/3' ix-y,p) s%p, {y,p). 



lim 



(4.25a) 



(4.25b) 



(4.25c) 



These equations are the most important result of this section. They simultaneously describe all "partonic" splittings, 
recombinations and scatterings from the distant past to the present. Note that these splittings and recombinations 
occur in all cells of coordinate-space. This is very different from the conventional approach where particles interact 
only when they are within ^/aroT of each other ||2^,^,0 . 



^'^A.k.a. the particle and anti-particle Wigner functions. 

^^The dg^ jx, p) in this section is the Wigner transform of the vector potential, y4,i^(a;,q), of Section^, as we demonstrate in 



Subsection [VE 
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Equations ( 4.25| ) are the phase-space QED analog of Mahklin's evolution equations A QCD version of the 



phase-space evolution equations should reduce to Makhlin's equations when integrating out the coordinate depen- 
dence. Geiger Q has derived a set of QCD transport equations based on Makhlin's work. While his derivation 
is very similar to our derivation of the phase-space evolution equation, he does use the gradient approximation to 



simplify his collision integrals. The QCD version of the transport equations we derive in Section IV F would reduce 
to his semiclassical equations if one applies this approximation. 



There are several ways to solve equation ( 4.25| ) but we propose only two methods in the following subsections. The 



first method is a perturbative scheme which we will use to derive the time-ordered version of the results of sections 
The second method is to derive transport equations from equation (4.25). 



E. Perturbative Solution to the Phase— Space Evolution Equations 



We can perform a perturbative expansion on equations (4.25) and get the leading contributions to the particle 

densities. We show this for both the photons and electrons surrounding a classical (scalar) point charge. 

We begin by stating the initial conditional for the particle densities and listing the other assumptions used here. 

— > > 

The initial electron and photon densities (at yo — > —oo) are s^^ {—oo,y,p) = dj:„ {—oo,y,p) = 0. We also take the 

initial scalar densities to be (— oo, y,p) — and 5< {yo, y,p) — M6{pq)S^{p — pi) S{p'^ — M'^)6^{ xop/ pQ — x). To get 
this form for , we localize the initial scalar wavepacket in momentum as discussed in Appendix |C 3[ In addition to 
assuming these densities, we must also neglect the mean field and drop the gradients in the scalar-photon coupling. 

1. Photons 

Since the scalar field only couples to the photons, the lowest order contribution to the photon den sity comes from 
the photons directly coupling to the initial scalar density. The cut diagram for this process is in Fig. For positive 
energy photons, the density is 

d<Ax,p) = [d^y^,^,D;^^,^,{x-y,p){27rfSHp-{qi+q2)) 
J (27r) (27rj 

xaemZ2(qi + g2)^'g< (y, qi) {qi + q2T' g< (y, 92) • 

Now, G^{x,p) — G^{x, —p) because ^ propagators obey the relation G-^{x,y) = G^{y,x). Thus, we can switch one 
of the g^ (y, q) to 5^ (y, —q), changing it from an initial state antiscalar to a final state scalar. Doing so, we have 

d^A^,P) = jd'y^,^,Dl^,Ax-y,p){2.)H\p^{q,-q,)) ^^^^^^ 
xaemZ^iqi - q2T'g< (y, qi) {qi ~ q2Y' 9^ (y, q2) ■ 

Now we can bring aemZ^{qi — 52)^ 9^ (y, qi) {qi — q2Y 9^ {Vj 92) into the form of the Wigncr transform of the scalar 
current. To do this, we take the final state scalar to be free and sum over all possible final momentum]^ Doing so. 



Eq. (4.26) becomes 

d^Ax,p)= I d^yD+ ,^,{x-y,p)J^i^^^^^Jy,p). 



Thus, dl]^ {x,p) can be identified with the A^i^(x,p) in equation (E 



Unlike Feynman perturbation theory, we can only specify the initial particle densities here. 
We perform this calculation in detail in Appendix DI 
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2. Electrons 



Since the electrons only couple to the photons, the lowest order contribution to the electron density com es fr om a 
photon splitting into electron-positron pairs. The cut diagram for this is shown in Fig. ^ From equation ( 4.25| ;) we 
have: 



'a/3 



y~ae,n{l'')a'a"S<„p„ {y , Qi) [Y ) 0" p' d"^^ (y, (72) • 
Using s<^ (x, q) = (a;, ~q), we find 



sif,^'0'i^ - y^p) (2^)'^'(p - (-91 + 92)) 



(4.27) 



Xaem(7^)a'a"s|-/„// {v , qi) il") 13" P' (^,92) ■ 



Taking the initial photon density from (4.26) and taking s^,,^,, {y, ^2) to be a final state positron, we recover equation 
(3.4). However here all of the propagators are retarded while the electron propagator in equation (3.4) is causal. 



F. Transport Equations 



In this section, we find a set of transport equations from the integral equations in (4.2£). We write the two equations 
of motion for the phase-space retarded propagator. Applying these equations to the phase-space evolution equations, 
we derive two sets of coupled integro-differential equations. The first set of equations are the transport equations and 
the second set are the "constraint" equations of Mrowczyhski and Heinz [PPOfI . The "constraint" equations describe 
the mass shift of the particles in medium. 

The equation of motion for the non-interacting retarded massless scalar propagator is 

d''G+[x) = 6^{x). 

The conjugate equation is 

d''{G+{x)y = 5^{x). 

Multiplying both sides of the first equation by (G+(?/))*, both sides of the second equation G^{y) and Wigner 
transforming in the relative space-time coordinate, we find two equations: 



{k + id/2fG+{x, k) = J d^x' e"'-'= {G+{x - x' 12))* 5'^{x + x' /2.) (4.28a) 
{k - id/2fG+{x, k) = [ d*x' e"'-'= {G+{x + x' /2)) S-^ix ~ x' /2) (4.28b) 



Inserting the retarded propagator in the energy-momentum representation (with mg — 0) and adding and subtracting 
the + and — equations, we find the equations of motion for the retarded propagator: 

2 

k ■ dG+ix, k) = -e{xo)S{x^)sm {2x ■ k) (4.29a) 

TT 

(aV4 - k^)G+{x, k) = -9{xo)S{x^)cos {2x ■ k). (4.29b) 

TT 

Taylor series expanding the sine or cosine and keeping only the lowest order is equivalent to performing the gradient 

expansion. 

Now, we apply the k ■ d and (5^/4 — k^) operators to the particle densities in equation ( f.25 ). On the right hand 
side, these differential operators act on the retarded propagators, so we can use their equations of motion to simplify 
the results. For scalars we get 
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P ■ dg^ {x,p) 



4 d'^qi d'^q2 2 ^ 
y 7^r~^ 7^r^ ~^(^o " 2/o)<5((x - y) ) 
(Z7rj (27rj ~ 



X sin (2(a: - y) ■ p){2nY5\p - (qi + 92)) 
Xa^mZ^qi +q2-id /2rg^ (y,qi) iqi+q2-id l'2Yd% {y^q^) 



(4.30a) 



+ d y -9{xo - 2/o)'5((x - y) ) sin (2(a; - y) ■ p)iQ^p{y) 

TT 

d^y -e{xo - ya)Si{x - yf) sin {2{x ~ y) ■ p)g^ {y,p) 

IT 



lim 



{d-'/4-k')g^{x,p)= I d'y 



4 d'^qi d'^q2 2 



~^(^o ^ 2/o)'5((a; - y)^) 
(27r) (27r) tt 



X cos {2{x - y) ■ p){27t)^6%p - {q, + 92)) 
xae™Z2(gi +q2~id /2)^5^ (y, 91) (gi + 92 - i a /2)''d|, (y, 92) 
+ / rf^'y -d{xo - yo)'5((a; - y)^) cos (2(x - y) • p)iQ^p{y) 



(4.30b) 



lim 



ya- 



d^y -e{xQ - yo)6{{x - y)^) cos (2(x - y) • p)g^ {y,p) 



Now because of the delta functions, the boundary conditions at yQ — s- —00 only contribute when \x — y| goes to 00, 
implying that we need g^ {x,p) as x — s- 00. The densities are zero here, so they drop out from these equations. 
The transport equations for the photons and electrons are 



p-dd%{x,p) 



TT 

d^qi d'^q2 ^o_\4x4 



d'^y 0{xo - yo)5{{x - yf) sin {2{x ~ y) ■ p) 



{2nf {2nf 



{2^)H\p - {q, + 92)) {ae^Tr [7^5^ (y, q,)-f''s^ (y, 92) 



(4.31a) 



+ aemZ^qi +q2 + id l2Yg^ (y, qi) {qi+q2 + id /2)"g^ (y, 92)} 



> 2 

X y" d-^y 6'(a;o - yo)'5((a; - y)^) sin (2(x ~ y) ■ p) 
d'^qi d'^q2 



{2^)H\p-{q^+q2)) 



{2nf {2nt 



(4.31b) 



X aem(7^)a'a"s|„^„ (y, ^l) il") f3" 13' d% (y, (72) 

+ (5„'/3' (±iS;|F(y))| . 

These equations almost have the form of the Boltzmann equation: the left side clearly is the Boltzmann transport 
operator and the right side is almost the collision integrals. If we were to expand the (co)sines in the collision 
integrals and keep only the lowest term, we would recover the collision integrals. Furthermore, if we were to do this 
same approximation to the QCD version of ( 4.31 ) we would arrive at Geiger's semiclassical QCD transport equations 

We also state the constraint equations: 



(974 - k')d% {x,p) ^ - j d% 9ixo - yo)S{{x - yf) cos (2(x - y) ■ p) 

7'"s^ {y,qi)'-f''s^ {y,q2) 
/2rg^ (y,(72)} 



in) (2n) 



d'^qi d'^q2 
(2^(2^ 

+ aemZ^iqi +q2 + id /2)''g^ (y, qi) {qi + q2 + i d 
+.9M^in|p(?/)| 



(4.32a) 
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X / d*y e{xo - yo)S{{x - yf) cos {2{x - y) ■ p) 

-^^(2vr)V(p-(,i+,2)) (4.32b) 
(271") (27r) 

X ae„r{Y)a'a"S^np„ (?/, ^l) iY)t3"p'd% (y, (72) 
+ (±»2^Mf(2/))} • 

The constraint equations give rise to the in-medium mass shift for the photons and electrons because, were we to 
derive the constraint equation for massive particles, we would find that (9^/4 — fc^) (9^/4 — k'^ + m?). Thus, the 
RHS of the constraint equations can be interpreted as an "in-medium" mass. Note that, even in the presence of this 
"in-medium" mass, particles scatter onto the light-cone. This is not a surprise since the particles are massless. Also 
note that the theta functions enforce the retarded time-ordering. Finally, we have not written the various constants 
in terms of there renormalized values. Dressing the particle densities by solving the evolution equations (which are 
nonpcrturbative) should be equivalent to using renormalized couplings. 



G. Implications for Parton Transport Theory 



The "source-propagator" picture must be valid for partons since the phase-space Generalized Fluctuation- 
Dissipation Theorem does not depend on the form of the self-energies. So, if we find the QCD self-energies and 
define the parton distributions appropriately, we could construct phase-space parton evolution equations. Of course, 
for these equations to have meaning, one must understand the role of Aqc d ■ ^qcd is usually interpreted as a mo- 
mentum cut-off in perturbation theory; as one approaches AqcD: nonperturbative effects increase and perturbation 
treatments break down. This interpretation may not be appropriate for several reasons. First, the phase-space evo- 
lution equations are non-perturbative objects, so there should be no cut-off in momentum. Second, it is not clear 
whether Age d should viewed as a cut-off in momentum or whether 1 / Aqcd should be viewed as a cut-off in coordi- 
nates. In fact, it may be that 1/Aqcd is simply a characteristic length scale for QCD bound states (i.e. hadrons) so 
that higher-order correlations (i.e. four and six point functions) need to be built into the Dyson-Schwinger equations. 
Presumably this would require some understanding of hadronization. Nevertheless, in the absence of a phase-space 
evolution equation, we can still use the Generalized Fluctuation-Dissipation Theorem as insight to build models. This 
is what we do in the next section. 



V. PARTON CLOUD OF A NUCLEON 



We cannot calculate the phase-space Parton Distribution Functions without a set of QCD phase-space evolution 
equations. Nevertheless, there is significant work calculating the Parton Distribution Functions in momentum-space 
and many of these results can be translated into phase-space. In particular, we show that the Leading Logarithm 
Approximation works in phase-space. Using the momentum ordering in the Leading Logarithm Approximation and 
a simple model of the nucleon we estimate the size of the sea parton distribution as a function of parton momentum. 



A. Parton Model and Leading Logarithm Approximation 



The parton model rests on a two simple assumptions: <C 1 (so perturbation theory is valid) and the parton 
lifetime is much larger than parton interaction time |2^. Both of these conditions are necessary to factorize a cross 
section in a momentum-space calculation. Typically the parton distribution functions are calculated using either 
DGLAP, BFKL, or GLR evolution, all of which are equivalent to applying a Leading Logarithm Approximation 
(LLA). In the LLA, we assume the parton is produced in a cascade represented by the ladder diagram in Fig. 
The probability of emitting the n*^ parton with longitudinal momentum fraction Xn and transverse momentum 
from this cascade is Q 

dP„ = ^^^. (5.1) 
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Thus, by ordering the momentum properly as we go down the ladder, we can pick up the largest logarithmic contri- 
butions to the n*'* parton's density. 

Most hadron colliders probe regions where the data are well described with parton distribution functions calculated 
within the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution scheme. DGLAP evolution is equivalent 
the Leading Logarithm Approximation in (LLA((3^)). New experiments at HERA are beginning to see evidence 
that Balitsky-Fadin-Kuraev-Lipatov (BFKL) type evolution is necessary to describe the parton distribution functions 
at small-a; BFKL-type physics is believed to be responsible for the rise in the number of partons as a; 0, 

however this rise can also be partially described by DGLAP-type physics BFKL evolution is equivalent 

the Leading Logarithm Approximation in 1/x (LLA(a;)). Unlike DGLAP and BFKL evolution, Gribov-Levin-Ryskin 
(GLR) type evolution does not have a simple momentum ordering because one sums terms with varying powers of 
1/x and l/q^ [^^. Because of the simplicity of the ladder structure and the momentum ordering needed to pick 
up the largest contributions, we will discuss both DGLAP and BFKL type partons in phase-space. 

We can apply the parton model and LLA in phase-space if both are modified appropriately. Assume that we are 
working in a regime where ^ 1 , so we can apply phase-space perturbation theory, and assume that all elementary 
particles are massless. Assume also that the probe is localized on the length scale of the parton cloud. This assumption 
is equivalent to saying the parton lifetime is large compared to the interaction time. 

Now, if we find the same singularities in both phase-space and momentum-space, then we know that the LLA 
will give the dominant contribution to the particle densities in phase-space. The Generalized Fluctuation Dissipation 
Theorem tells us that the parton density has the form 



7^(x,p)= I d^yG+{x-y,p)i:^{y,p). (5.2) 

th 



The self-energy, I]<, is given by the parton ladder in Fig. hB and the n segment of I]< is shown in Fig. hM In 



momentum-space, the cut rung gives a d'^k/\ko\ which leads to the dx/x in equation (5. 1^ To see how the factor of 



d'^k/\ko\ arises in phase-space, one needs only look at the electron source in Section [IIB 1. The electron source has 
exactly the form of the segment in Fig. ^ and in that calculation we found exactly this factor of d'^k/\ko\. The fact 
that we find the same factor of d'^k/\ko\ in both the energy-momentum representation and in phase-space simply 
reflects the fact that the cut parton density is 9{kQ)5(k^) in both cases and we sum over final parton states. The 
factor of dq^ /q^ occurs in equation (|5.lD because of the integration over the leg's propagator \/q^. In phase-space. 



the 1/q pole is tied up in the Wigner transform of the retarded propagator, but it is still there: 
G^{x,qj = / 7—4 e 



{2nf (q + q'/2r + ^<qo + 9^/2) {q - ^ jlf - ie{q^ - q'j2) ' 

Thus, this segment of the parton ladder produces the same divergencies in phase-space and momentum-space. What- 
ever orderings are needed to produce the leading contributions in momentum space will produce the same leading 
contributions in phase-space. 

Our self-energy has the same ladder structure as the electron source in Section |l|, so we know the spatial structure 
of the n'^ parton's source is given by the n — I*'' parton's distribution. Iterating back to the 0*'' parton (a valence 
quark), we see that the shape of the valence distribution sets the shape of the sea parton source. We take the 
valence quark wavefunction to be uniformly spread throughout a bag with radius Rbag- Since we are interested in 
high-energy collisions, we take the nucleon bag to be moving to the right with 4-momentum = (Pqi Plt^t) with 
Po ~ Pl ^ Mn. Thus, this nucleon has 4-velocity {1,vl,0t) and the bag is contracted in the longitudinal 

direction by a factor of 7 = 1/^1 — ^ 1. We assume the partons lose memory of the original valence quark 
momentum as one goes down the ladder. Thus, any momentum/coordinate correlations in the source function should 
be washed out by the spatial integrations in equation ( |5.2| ). One might expect that the sea partons forget the shape 
of the nucleon bag as well, but we show that the partons cannot propagate far enough from the original source for 
this to happen. 
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B. Large-Q^ (DGLAP) Partons 



In the large-Q^ regime, the parton density is low but as{Q^)\a{Q^ / K^(jjj) > 1. Here the largest contribution to 
the leading log ladder comes from large logarithms.]^ To get the largest contributions from these logs, we order 
the momenta as we move down the ladder: 

-ql » -ql-l » • • • » -ql » l/Rlag ~ ^QCD- 

Here qf is the virtuality of the i*'* leg. The kinematics at each leg-rung vertex ensure that the momentum fraction 
carried by each leg is also ordered: 

1 > Xi > . . . > Xn^i > Xn- 

Whether a rung or leg is a quark or gluon is irrelevant, provided fc^ = and the ordering holds. Now, given that 
the proton has longitudinal momentum and the rungs and legs are massless, each generation of partons must have 
energy qno ~ XnPh and transverse momentum of q^ w —q^ <C x^P^^. 

Let us discuss the general features of the parton cloud. The retarded propagator lets the n*'* parton propagate out 
to Rn± ~ he/ y^—q^ transverse to the parton momentum and to i?„|| ^ hc/qno — hc/xnPh parallel to the the parton 
momentum. Since xPl ^ pT, the parton momentum is approximately parallel to the nucleon momentum. Because 
Rni_ ^ Rbag, the partons can never get far from the bag in the transverse direction and transverse spread of the 
partons will dominated by the bag size: Ai?T ~ Rbag- On the other hand, the longitudinal spread of the partons 
is roughly given by Ai?i ~ Rbag/ J + fic/xPL, so can be dominated by the longitudinal propagation distance if 
X <C MpfRbag/fic. In fact, for very small x (i.e. x ~ AlNRbag/jhc) the spread of the partons can meet or exceed the 
nucleon bag radius. Furthermore, the actual distribution may be somewhat broader due to the propagation of the 
virtual partons between the subsequent emissions along the ladder. 

So, in our picture, the sea quark and gluon large-Q^ distributions have the same transverse size as the parent 
nucleon, but the longitudinal size can be significantly bigger than the parent. Furthermore, the drop off in the parton 
density in the longitudinal direction occurs at the characteristic radius of ~ Hc/xPl. This picture of the nucleon is 
consistent with the uncertainty principle based arguments of A. H. Mueller ]20|] , later user by Geiger to initialize the 
parton distributions in his Parton Cascade Model . 

C. Small-s (BFKL) Partons 

In the small-x regime, the parton density is high and as{Q^) ln(l/a;) > 1. The small-x partons are mostly gluons. 
In this regime, the leading logs come from the l/cc-type singularities, i.e. from the cut rungs. Typically the rungs 
are thought of as "reggeized," meaning that each rung in the ladder is allowed to radiate gluons, leading to a large 
multiplication in the parton density. Since leading logs come from the l/x type singularities, the largest contributions 
come about by strongly ordering the longitudinal momentum fraction as one moves down the ladder p^ : 

1 » a;i > . . . > Xn~l > Xn- 

BFKL-type evolution has only a weak dependence on the virtuality of the partons as we move down the ladder, so 
we assume to be fixed: q^_i ^ qn^ ^/Rbag- This does not significantly effect the results of the analysis [|4|. 

Now we must understand how the transverse momentum and energy of each parton leg changes as we go down 
the ladder. A well known effect of iterating the BFKL kernel (equivalent to moving down the ladder) is that the 
transverse momentum undergoes a random walk in ln{q^) H^Js^. In fact, after iterating through a sufficiently large 
number of rungs, the spread in the qT distribution is given by: 



^^Q^ can be taken as the typical momentum scale of the process. In the case of a DIS probe, this is the momentum transferred 
by the probe. 
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where C = ^^28C(3) = 32.14a^. Thus, qlj. can be orders of magnitude larger or smaller than qfrj^. We restate this 
as 

qlr-qlre^'-'^^^^^^- (5.3) 

We will consider the extreme cases of the transverse momentum and comment on the typical case, q'^rp ^ q^rp. 

If the random walk results in a large transverse momentum, we will have qf^j, ^ qfj, ^ — ~ [xuPlY ■ Thus, the n*-^ 
parton will have 3-momentum in the transverse direction. We know that the parton can only propagate to a distance 
of roughly R\\ ^ /ic/|go| = fic/ \/ q^ + q^ + {xPlY the direction parallel to q. Since he/ yjq^ + q^ + {xPi^Y w 
hcl\qT\ ^ Rbagi the parton cannot travel far from the original source in the transverse direction. On the other 
hand, the parton's longitudinal spread can be larger than the longitudinal bag size. The parton can propagate to a 
distance of R± ~ hcj \J —(p in the direction perpendicular to g, so we can expect a longitudinal spread of the parton 
distribution of Ai?^ ~ Rbag/l + ^c/ ^/—q^- Since Rbag ^ ^c/ \/—q^, this additional spread can not match the spread 
of the DGLAP partons. 

If the random walk results in a small transverse momentum, we will have <^ qfrp ^ —q^ ^ {xuPlY'- In this 
case, the n*'* parton will have 3-momentum in the longitudinal direction. As in the case of the DGLAP partons the 
additional transverse spread is Ai?^ ^ fi.c/-\/— g^ <^ Rbag and so is negligible. The additional longitudinal spread is 
Ai?i ^ ?ic/|go| — he/ ^Jq^ + {xPlY + may be significantly larger than the bag radius because the parton is 

space-like. 

Summarizing both possibilities, the BFKL parton distributions have the same transverse spread, Ai?^ ~ Rbag, 
but different longitudinal spreads. The longitudinal spread may range from ARl ^ Rbag/l + fi.c/-\/— g^ <C Rbag to 
ARl ^ Rbag/l + he/ \/ q^ + [xPlY + 9t ^ Rbag for partons with space-like momentum. The fact that the spatial 
extent of the BFKL cloud is so large in the longitudinal direction suggests that the small-x partons (which are mostly 
gluons) can see the color charge of any other nucleon in the longitudinal tube centered on the parent nucleon. This 
suggests that we should treat the nucleus as a whole as a source of color charge in the spirit of McLerran-Venugopalan 
model ||]. 

The large longitudinal extant of the small-x cloud has another consequence: in a zero impact parameter nucleon- 
nucleon collision, we would find that the soft (BFKL) partons interact much earlier than the harder (DGLAP) partons 
because of their greater longitudinal spread. This, coupled with the large density of small-a; partons, leads to earlier 
entropy production and stopping of the soft partons. 

VI. CONCLUSION 

We have made progress toward specifying the initial phase-space parton distributions of a relativistic nuclear 
collision. Regardless of the kinematical regime, the transverse spread of a parton distribution is dominated by the bag 
radius ^ 1 fni. The longitudinal spread of a parton distribution varies from roughly ^ Rbag/l + he/xPL for moderate 
to large x (i.e. for DGLAP partons) and from ARl ^ Rhag/l+'he/ \J —(p to ARt ^ Rbag/l^-'he/ \/q^ + {xPlY + g|, 
for small x (i.e. BFKL partons). Since the small x partons have a large longitudinal spread and a high density, we 
expect the small x partons to interact much earlier than the large x partons in a typical nuclear collision. This would 
cause earlier entropy production and higher stopping than one expects in models that include only DGLAP parton 
distributions such as the PCM HIJING |], and others. 

Using time-ordered non-equilibrium methods, we derived phase-space evolution equations for QED, illustrating 
how to find them in QCD. Unlike conventional transport approaches, our calculation does not rely on the gradient 
approximation. Thus, it should work on all length and momentum scales. These phase-space evolution equations 
describe the complete evolution of a system from some time in the distant past to the present, including all splittings, 
recombinations and scatterings of the particles. One can use these evolution equations perturbatively or to derive 
semiclassical transport equations. These evolution equations rely on the Generalized Fluctuation-Dissipation Theo- 
rem. This theorem states that a particle's density is the convolution of the Wigner transform of its self-energy and 
a phase-space propagator. The Generalized Fluctuation-Dissipation Theorem is quite general and can be directly 
applied to QCD. 

In conventional Feynman perturbation theory, we found the reaction rates (and hence the cross sections) can be 
written in a parton model form. In other words, they take the form of a reaction rate density convoluted with a 
phase-space Parton Distribution Function. This phase-space PDF is the parton number density and has the form of 
a phase-space source folded with a phase-space propagator. Our work with the Weizsacker-Williams Approximation 
demonstrates that the Parton Distribution Functions can be defined in phase-space. 
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In order to illustrate how the propagators and sources work in phase-space, we calculated the effective photon and 
electron distributions. We found that both the retarded and Feynman propagators propagate pa rticl es to distances of 
= hc/mhi{\qo\, \q\) parallel to the particle's momentum and to distances of ^ R± — hcj perpendicular to 
the particle's momentum when ^ 0. When <^ = 0, the particles tend to follow their classical paths with deviations 
from this path being of order l/|go|- Furthermore, the retarded propagator can only send particles forward in time 
and inside the light-cone while the Feynman propagator sends particles both forwards and backwards in time and 
both inside and outside of the light-cone. We also described a phase-space source that included a simple "partonic" 
splitting: the electron distribution of a point charge. These electrons are created when a virtual photon splits into an 
electron-positron pair; the diagram for this process is the first segment of a parton ladder. Wc found that shape of 
the electron's source is controlled by the parent photon's distribution. 

We hope that we have provided insight into the behavior and calculation of the phase-space densities. Specifically, we 
hope the "source-propagator" picture of the Generalized Fluctuation-Dissipation Theorem and the resulting phase- 
space evolution equations can be coupled with appropriately defined phase-space parton densities. The resulting 
theory could describe the various many-particle effects we expect in a nuclear collision at RHIC or the LHC and it 
could incorporate parton model phenomenology. 
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APPENDIX A: ELECTRODYNAMICS WITH FERMIONS AND SCALARS 



Throughout this paper we use QED and Scalar Electrodynamics to describe the interactions between the electrons, 
photons and scalars. In this appendix, we review the QED lagrangian, equations of motion, equal time commutation 
relations and contour Feynman rules. The lagrangians, etc., for QED and for Scalar Electrodynamics are given in 
many places [^^. Even so, we restate them here both to keep this work self-contained and to clarify our notation. 
We do not include the renormalization counterterms nor the gauge fixing terms for the photons although they can be 
easily included. We work in the Lorentz gauge. 

The lagrangian for scalar QED coupled with spinor QED is 



+ {df,(j)*ix)){df'c^{x)) ~ M^*{xmx) 



(Al) 



-eip{x) 4-{x)i;{x) - iZeAt'{x){(j)*{x) ^x)) + Z^aemA^{x)(t)* (x)(t){x). 

Here ipaix) is the fermion field, (f>{x) is the complex scalar field, A^(x) is the photon field, and F^^(x) 
d,^Afj_{x). The masses of the fermion and scalar fields are nie and M respectively. The electrons couple to the photons 
with strength e while the scalars couple with Ze. 

The second quantized field operators satisfy the standard equal time commutation relations: 



d^,At,{x) 



Ap,{t, x), x ) 
Ap,{t, x), Ap^(t, x ) 

■4;a{t,x),1pj;{t,x') 

${t,x),${t,x') 
(t){t,x),^{t,x') 



4TTg^^6^{x- x') 
Afj,(t, x) , Ai^(t, X ) 

OaffS^ix- x') 

S^{x-x') 
4){t,x),(p{t,x') 



= 



= 



(A2a) 
(A2b) 
(A2c) 
(A2d) 
{A2e) 



The lagrangian ([Al[) leads to the following equations of motion: 
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infix) = d^F'^'^ix) 

= (i ^ — e ^{x))'4'{x) — meipix) 

= (5^ + ieA^{x)){d'' + ieA" {x))(p{x) + M% 



(A3a) 
(A3b) 
(A3c) 



The electromagnetic current operator is 

= e%l}{x)^^ip{x) + iZe(l)*(x) 



x). 



These equations are solved by the Green's functions in Section [VA in the limit as e ^ 0. 
We list the contour Feynman rules for spinor and scalar QED below: 

1. The vertex Feynman rules are summarized in Table |. 

2. The contour propagators are summarized in Table ||. 

3. Every closed fermion loop yields a factor of (—1). 

4. Every single particle line that forms a closed loop or is linked by the same interaction line yields a factor of iG"^ . 
Notice that the second scalar coupling has higher o rder t han the rest of the couplings. So we neglect this coupling in 



the derivation of the evolution equations of Section IV D 



APPENDIX B: THE CROSS SECTION IN TERMS OF PHASE-SPACE DENSITIES 

In this appendix, we discuss the cross-section in terms of phase-space quantities. Since the cross section is measured 
by scattering a beam of particles off a target, we define the cross section in terms of the projectile/target reaction 
rate density and the projectile flux. The beam is uniform in the beam direction and in time on the scale of the 
projectile/target interaction. Thus, the beam can only directly probe the transverse structure of the interaction 
region. Even this transverse information is washed out in the typical experiment, since the beam is usually uniform 
in the transverse direction on the length scale of the interaction. In the limit of a transversely uniform beam, we 
recover the conventional definition of the cross section. Since we consider only simple scattering problems, we work 
in Feynman perturbation theory where we can specify both the initial and final states of the reactions. 

The beam is a collection of single particle wavepackets distributed uniformly throughout the transverse area A of 
the beam. For the sake of illustration, we take these particles to be scalars. The Wigner function of these incident 
wavepackets is 



f{x,p) 



1 



2VpqJ {2TrY 



e-^-P f{p + p'/2)r{p-p'/2) 



where the wavef unction f{p) is given hy 



d'^p 



fip) \P) ■ 



(Bl) 



(B2) 



We will assume the beam to be uniform in the longitudinal direction with length L and to be turned on for macroscopic 
time T. The quantities A, T, and L are much larger than the projectile/target interaction region. 

The projectile/target interaction region is characterized by a reaction rate density Wi^/(x). We assume the reaction 
rate density to be localized in both space and time. This reflects the small spatial extent of the target and the short 
interaction time compared to the beam lifetime. The reaction rate is trivially related to the reaction probability: 



d^xW.^fix). 



Thus, the reaction rate is easily identifiable in the calculations in sections |l,III. For example, in the process 7B 
in Fig. |l|b, the reaction rate density is WjB^B'ix, q). For the process AB A'B' in Fig. [l^., it is 



(B3) 



B' 



The delta function that puts the particle on-shell is absorbed into f{p) 
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Wab^A'B'{x) = / d^r-^J'X'{x + r/2)Dl,^,Ar.q)J'B'' {^-r/2). (B4) 



Note that the reaction rate density is a function of the average space-time location of all the vertices in the process. 

The cross section is the effective area of the target, so we define the cross section as the integral over the beam face 
of the fraction of incident particles that interact with the target per unit area: 

[ cP f ^ scattered particles \ / ("^ incident particles \ 
J \ unit area / / V unit area / 

The number of incident particles per unit area crossing the target plane is the particle flux: 



# incident particles 
unit area 



L/2 

= J^inc I dxL fi ■ j{x) = T{xt). (B6) 

L/2 



Here n is a unit normal to the target plane and Mine is the number of particles in the beam. The single particle 
current is given in terms of the incident particle Wigner function by 

J(f) = j £pdp^vf{x,p). (B7) 

We need not average over time because the beam is uniform on the time scale of the reaction. The number of scattered 
particles per unit area is found by multiplying the number of incident particles by the reaction probability per unit 
area: 



# scattered particles 



L/2 i.T/2 



M,nc dxL dxoW,^f{x)=J\f,ncm^f{xT). (B8) 

unit area J-l/2 J-t/2 



Thus, the cross section is 



,2 M«c Wi^fixr) 



d'xT '"''-7i'^' . (B9) 

J^(Xt) 



In equation (B9), all longitudinal and temporal structure of the interaction is washed out by the beam. Furthermore, 
in a any practical experiment, the wavepackets are delocalized in the transverse direction on the length scale of the 
interaction region. Thus, the transverse structure of J-{xt) is gone and the flux reduces to T = A/i„cKI/^: where \v\ 
is the mean projectile velocity. The flux can then be pulled out of the transverse integral in (^9|). The transverse 
integral of the reaction probability per unit area is A/'mc|<5'i_+/p, so the cross section becomes 

This is the conventional momentum space cross section in our choice of normalization. 

APPENDIX C: WAVEPACKETS 

Throughout this paper, we use wavepackets in the initial and final states of a reaction to provide spatial localization 
or delocalization. In this appendix, we detail the construction of an initial or final state wavepacket and discuss the 
limits of cither a completely localized or delocalized wavepacket. 

1. On— Shell Gaussian Wavepacket 

An initial (or final) state ket can be written with wavepackets: 

-^f{p)\p)- (CI) 
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The corresponding Wigner function of the particles is 



(271-) 

2VpoJ {27rf'' 



f{p+p'/2)r{p-p'/2). 



(C2) 



Particles in either the initial or final states are on-shell, so they can be expanded in momentum eigenstates. We 
choose our wavepacket to be a Gaussian superposition of momentum eigenstates with a momentum spread a: 

cl>{p) = N5{p^ - M2) exp [-(p - Pif/2a^] 

The Wigner transform of this wavepacket can not be done analytically except in the limit when \pi | » a. In this 
limit, Pi Ri p so our wavepacket is localized in momentum giving the following Wigner density of particles: 



f{x,p) 



87rp2 



(5(p2-M2)exp 



■Pi? 



{2<jy/2^f exp [-2(t'^{vxo - xf] . 



(C3) 



Here v = p/po is the velocity of the wavepacket. Thus, the particle's Wigner function is a Gaussian in both momentum 
and space. The spread in momentum is the inverse spread in space. The centroid of the Gaussian follows the particle's 
classical trajectory. The energy of the packet is set by the delta function out front. We have not constrained the 
particle in energy so this density contains both positive and negative energy contributions. 



2. Delocalizing the Wavepacket in Space: Free Wavepacket 



In accordance with the uncertainty principle, the wavepacket becomes completely delocalized in space in the limit 
of complete localization in momentum (i.e. a ~> 0). In this limit, the spatial Gaussian approaches unity and the 
momentum Gaussian becomes a delta function. After working out the normalization, we find 

zvpo 

This is no surprise since we squeezed the state into a momentum eigenstate. 



3. Localizing the Wavepacket in Space: Classical Wavepacket 



A classical particle is localized in both space and momentum, a seeming violation of the uncertainty principle. In 
real life, this is not a problem since the reason classical particles appear localized is that we probe them on length (or 
momentum) scales too coarse to rcisolvc; the; interesting qiiantum features. In the case of our Gaussian wavepacket, 
this amounts to probing the distribution on length scales much larger than l/cr. In this case, the space Gaussian is 
too localized to resolve and we can replace it with a delta function. Additionally, we assume that a is large, so we 
can replace the momentiim Gaussian with a delta function as well. 

Making these approximations, we find the Wigner density of a classical particle: 

/-'™i(ar,p) = ^{2wfS^{p- Pi)5{p^ - M^) 5\vxo - x) (C5) 

Here we have inserted the correct normalization for the wavepacket. This density corresponds to an on-mass-shell 
particle that follows its classical trajectory vxq = x. Again, we left in both positive and negative energy contributions. 



APPENDIX D: THE CLASSICAL CURRENT 



In this appendix, we derive the classical current used in the effective photon distribution calculation. For the sake 

of illustration, we take our point particle to be a scalar particle. The derivation goes in three steps: first we define the 
Wigner current of a scalar particle, then we derive the photon/scalar interaction vertex in phase-space, and finally 
we localize the initial and final states of the scalar to give the classical current. 
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1. Wigner Current 



We begin by restating equation (2.3): 



^r(^,9)- / 4^,e-^^-^{A'\f{q + q/2)\A){A\j^^q-q/2)\A'). 



(Dl) 



We write the initial and final state bra's and ket's according to equation (B2). Rewriting equation (Dl) in terms of 
initial and final Wigner densities, 



J'^{x,q) 



77717 77r74 {x,Pi) {x,pf) {2tt)'^S^{p, - Pf - q) T ^^{q,p.,,pf). 
(27r) (27r) 



(D2) 



We assume that the initial and final wavepackets are localized in momentum and some-what delocalized in space. 
Shortly, we assume that we probe the current on length scales much larger than even this delocalized space distribution. 



2. Scalar Vertex 



T^^{q,Pi,Pf) is not quite th e W igner transform of the ^AA' vertex, although it does arise from performing the 
Wigner transform in equation (Dl). It is defined by 

{27:rS\p,~Pf - q)r^.(q,p„p/) - W'pjp^ f {pf\j^{q + q/2) \p,) {p.,\jUq - q/2) \pf) (D3) 

J (27rj 

Using the matrix element 

{Pfljt^iQ) \Pi) = eZ{2TT)^S'^{p^ ~pf ~q) 



iPfP, 
we get 

^t^A^^Pt^Pf) = a^rnZ'^iPi+Pf + ]^{Pl+Pf))^,{p^+Pf - l^{Pi+Pf))u (D4) 

The relative momenta, pi and p/, become derivatives on x in the current (p2|). We assume the wavepackets to be 
uniform the the reaction's length scales, so we ignore the derivatives and arrive at the phase-space scalar vertex 

^tiv{,q,Pi,Pf) = Olem.Z'^{p^+Pf)^{pi +Pf)u- (D5) 



3. Classical Current 



We are now in a position to derive equation (2.19) for the classical current density in phase-space. First, we take 
the final state to be a momentum eigenstate and sum over it. Since the final state is localized in momentum around 
Pf, this is not a bad approximation. Second, we take the initial state to be a classical wavepacket. In other words, we 
assume that the initial state is localized in momentum and delocalized in space but we probe it on such large length 
scales that we still see a spatially localized wavepacket. So, putting equations (C4), ( |C5| ) and (D5) into (D2) and 
summing over final states, we get 

J^''{x,q) = 2'KaemZ'^ v^v^ 6^{x - Xov)pio5{{pf + qf - Af^). 

Using p^j: = and « Pffj./Pio and assuming q^/pio ^ g • we get the classical current: 



J 



:lassical 



(D6) 



Note that this current allows for emission of both positive and negative energy photons. To use the retarded propa- 
gators in Section 1^, we need a 0{qo) in equation (DC). We can do this by suitable choosing pf and pi and restricting 
the initial and final states to have only positive energy. 
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APPENDIX E: PHASE-SPACE EFFECTIVE PHOTON DISTRIBUTION OF A STATIONARY POINT 

CHARGE 



In this appendix, we describe limit of u = of the photon distribution of the point charge in Subsection fl] 
Since the spatial dependence of the effective photon distribution is controlled by the Wigner transform of the vector 
potential, Afj_{x), we only discuss Afj_^{x,q) here. When v — 0, the photon vector potential becomes Afj_{x) — (|f|,0) 
so A^^{x, q) is the Wigner transform of the Coulomb potential. 

We take the point charge to be resting at the origin and emiting photons with four-momentum g^, = (go, q}- Putting 
w = in Eq. ( |2!27| ), wc find 



1 



^oo(a:^,'?) =32TT'aemS{qo)Y^Ai2\xMcos{0),2\x\\q\sm{9)) 



A,, 







191 



(El) 



where 9 is the angle between x and q and the dimensionless function A is given in Eq. ( ^.28D . Clearly the photon field 
is time independent and is composed of only zero energy photons. Furthermore, by virtue of the l/\q\ singularity, the 
photon field is mostly composed of low momentum photons. 

In Fig ^ we plot the dimensionless function ^ as a function of x in the plane defined by x and q. Note that 
the central region of the distribution is circular, but becomes elliptical as one moves away from the center. In the 
transverse direction (i.e. the direction perpendicular to the photon three-momentum), the distribution approaches 
zero, but never goes negative. The width in the transverse direction is approximately 250 fm. In the longitudinal 
direction, the distribution drops to zero at about xl « 250 fm and oscillates about zero for larger distances. These 
oscillations are expected for a Wigner transformed quantity and simply reflect the fact that xl and qL are Fourier 
conjugate variables. 

Because the photon source is a point source, th e the shape of the Coulomb distribution comes directly from the 
shape of the the retarded propagator in Subsection II D 2. Thus, we ca n estimate the width of the photon distribution 
using the estimates of the retarded propagator in Subsection II D 2. In the both the longitudinal and transverse 
directions, the propagator width is ^ hc/\qL\ — 250 fm, which is exactly the width we measure from the plots. 



APPENDIX F: EFFECTIVE ELECTRON PHASE-SPACE DISTRIBUTION WITH Me / 

In this appendix, we calculate the effective electron distribution for electrons with a mass much larger than their 
momentum. This calculation is not included in Section III because it is not relevant for partons. 



When the electron momentum is much smaller than its mass, we can use Remler's propagators for massive particle 
p3[ . His propagator takes one of two forms depending on whether the electron momentum is space-like or time-like. 
His propagator is discussed in Appendix We show a sample electron density for both the time-like momentum 
and space-like momentum cases. The momenta of the electrons are chosen to satisfy the requirements that p ■ v > 
and fcTmax be real. These requirements are equivalent to the requirement that p • v > me/7. Since p <C rue, we must 
also have 7 ::g> 1. 



In Section [II B 1 we show that the electron's source is controlled by the parent photon distribution so we show the 



parent photon distributions next to the electron distributions in all subsequent plots. We do not restrain the photons 
to have go > as in Section III, so our sources include negative energy contributions. Because we use retarded photon 
propagators in our source, these calculations only serve to illustrate how Remler's propagators function. In fact, 
had we restricted qo > 0, there would not be enough momentum-space to perform the kx integrals and the electron 
distribution would be zero. 



1. Feynman Propagator for Particles with Space— Like Momentum 

The propagator for electrons with space-like momentum is: 

P \ 1^1 1 



G%Ax,p) 



dr Ax - 



-2me|r| _ 



p2 J 2'mey'—p^(m'^ — p'^) 

|-\/— p2 COS (2r-\/— p2) 4- sin (2|r| -v^— . 



38 



Remler's propagator for space-like electrons has a very simple interpretation. First, the delta function forces the 
electron to follow its classical trajectory, but with the electron velocity defined as = Pfj.ly—p^- The exponential 
in proper time strongly damps propagation that extends farther in time than l/2me along the classical trajectory. 
The fact that the proper time can extend forward or backwards in time simply reflects the boundary conditions of 
the Feynman propagator. Next, the sine and cosine cause the expected Wigner oscillations. The rest of the terms 
simply give normalization. Finally, this propagator allows propagation outside of the light-cone, but such propagation 
is strongly damped. This may seem strange, but should come as no surprise: the coordinate space propagator for 
massive particles propagates particles outside of the light-cone [^ . 



2. Massive Electrons Distribution for Electrons with Space— Like Momentum 



We now perform the integrals over (PkT and d^x in equation (|3.4|). The d^x is a trivial delta function integral and 
the integral over (PUt can be done numerically. On the left in Fig. |18|we have a sample cut through the phase-space 
density for electrons with a typical space-like 4-momentum {p^ = (0.05, 0.008, 0.06, 0.0) MeV/c). On the right is one 
of the underlying photon distributions. We chose the source velocity so that 7 = 12.47. This velocity is a compromise 
between having enough momentum space available for the electron and rendering the plot unreadable because of the 
Lorentz contraction. 

Now we examine these plots. First, we see the contributions from retarded emission and propagation (upper left 
electron pancake) and from advanced propagation (lower right pancake). Let us concentrate on the retarded electrons. 
At some time in the past, the photons split into the electrons and the positrons. The electrons than propagate forward 
along their classical trajectory until they reach the location of the left pancake. Notice that this pancake has nearly 
the same size as the photon distribution on the right. The other photon distribution, corresponding to the other root 
of positron momentum, has a slightly different tilt and width, but the difference in the plots is not noticeable. The 
electron pancakes are slightly larger than the photon pancake, presumably because of momentum broadening from 
the emitted positron. The advanced electrons have exactly the same shape and size as their retarded brethren, but 
they followed a time-reversed classical trajectory, coming from some time in the future. 



3. Feynman Propagator for Particles with Time— Like Momentum 

Remler's propagator for massive particles with time-like momentum is 



1 



2rne vP^ 



sin(2T(Vj3^-me))^4/^_^_ p 



{\/p^ - me) 



_ sin(2r(vV +me)) ^4 I ^ _p_ 

{y/p^+nie) \ \fl 

Remler's time-like propagator does not have as simple an interpretation as his space-like propagator. Delta functions 
still keep the particle on its classical trajectory, but the integrals in "proper time" are Fourier sine transformed along 
this classical trajectory. Thus, a simple peak in the underlying photon distribution will get Fourier transformed into a 
series of peaks and valleys in the electron distribution. Furthermore, the advanced and retarded branches enter with 
different signs, so we have large negative contributions from the advanced branch. 



4. Massive Electrons Distribution for Electrons with Time— Like Momentum 

Despite the difhculty in interpreting the propagator, the d^x and d^fc^ integrals can be done. A sample cut 
through the phase-space distribution is shown in Fig. These electrons have a typical time-like 4-momentum 
(pp = (0.05,0.005,0.04,0) MeV/c) and the source has a 7 = 12.47 and is moving to the right. Again, the source 
velocity was picked as a compromise between readability of the plot and available momentum space for the electron. 

We see the difficulty in interpreting the electron distribution, the Fourier transform took a simple photon peak and 
produced a series of large electron peaks. The retarded branch corresponds to the envelope of large positive peaks 
on the upper left. The advanced branch corresponds to the envelope of large negative peaks in the lower right. Each 
peak in the pair of envelopes appears to be a Lorentz pancake, but the envelope as a whole is significantly broader 
than the underlying photon peak. Presumably, averaging the distribution over unit areas in phase-space would result 
in a much tighter distribution. 
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APPENDIX G: FREE SCALAR PROPAGATORS IN PHASE-SPACE 



In this section, we state all of our phase-space propagators, discuss the symmetries of the massless propagators 
in some detail and outline the derivation of the retarded and Feynman scalar propagators. The massive Feynman 



propagator is discussed by Remler |13| so our discussion is brief. The Dirac and vector propagators differ from the 
scalar propagators by the inclusion of either spin projectors (in the case of Dirac particles) or polarization projectors 
(in the case of vector particles) so we do not need to discuss them. 
We define the Wigner transform of any translationally invariant propagator as 

G{x,p) = [ ^^e~^--P'G{p + p'/2)G\p~p'/2) 

J (27r) (Qi) 



d'^x' e"'-fG(a; + x' /2)G\x - x' 12) 



The vacuum propagators that we use are 



{p) ^-{p^~m^± iepo) ^ (G2a) 
(p) ^-{p^-m^± le) " ' (G2b) 



1. Massless Scalar Propagators 

a. Symmetries 

A time reversal transform in coordinate space is equivalent to a simultaneous reflection in time and energy in phase- 
space. Under time reversal the + and — propagators change into one another while the Feynman and anti-Feynman 
propagators remain unchanged: 

G+(xo,f,po,p) = G~(-a;o,f, -po,p) (G3a) 
G''{xq,x,pq,p) = G''{~xo,x,-po,p). (G3b) 

A parity transform in coordinate space is equivalent to a simultaneous reflection in a space coordinate and the 
corresponding momentum coordinate. Under a parity transformation, all of the propagators remain unchanged: 

G=^ {xo ,x,Pq,p) = G^ {xq ,-x,Po,~p) (G3c) 
G''{xo,x,pq,p) = G''{xo,-x,po,-p). (G3d) 

The Feynman propagators have another (rather amusing) argument switching symmetry. Here all the space-time 
components are switched with the the corresponding momentum-energy components: 

G^(a;,p) = G"(p,a;) (G3e) 

Finally, the Feynman and anti-Feynman propagator are related through a complete reflection of all of the space or 
momentum coordinates: 

G%x,p) ^G-{-x,p) (G3f) 
G'{x,p) ^G-{x,-p). (G3g) 
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b. Propagators 



We now present the massless Feynman and retarded propagators. The advanced and anti-Feynman propagators 
can be recovered using the symmetry relations above. Since all of the massless scalar propagators are dimensionless 
and Lorentz invariant, we expect that they will be functions of x ■ p and x^p'^ and possibly theta functions in energy 
or time. In fact, the propagators are: 



[ sgn(a;^) + sgn(p^) + 2 sgn(x • p) 



1 

47r 



G'ix,p) 



G+{x,p) ^lg(xo)g(x^)g(A^) '^"^^^^ 



(G4a) 



(G4b) 



Here the Lorentz invariant A'^ is given by A^ = (x • p)^ — x^p^. Since we discuss how the propagators work in the 
Sections [ID 2 and [II B 3, we do not do so here. 



c. Derivation of 

The Wigner transform of G+ is easiest to do in coordinate space. In coordinate space, G+(x) — j^9{xo)S{x^), so 
t he Wigner transform integral is a series of delta function integrals. Performing the first delta function integral in 
(Gl), and simplifying the theta functions, we find 



G+{x,p) = 



0{xo) 



° dx'Q^ Ax"^ + x'^^ I d!^s'e"'-f(5(a;' • x). 



Using 2'k5{x) = dae"", we can do the angular integral, giving us a Bessel function: 



G+{x,p) ^ AT:e{xo)e{x^) / dae'°"^Ja{iyJl-a^) 



Here -q = 2{pi^\x\ — zoi • p) and ^ = 2\/ x '^jfp- — (p- xW). This integral is in any standard integral table [Q. After a 
bit of simplification, one gets the result (G4b). This result can be checked by performing the Wigner transforms in 
momentum space, but the contour integrals needed for this calculation are quite tedious. 



d. Derivation of 



The simplest derivation of G^{x,p) is far more complicated than the derivation of G^{x,p). We start by finding 
the transport-like equation of motion for the Wigner propagator.]^ The derivation is simple and very similar to the 
derivation for the retarded equation of motion in Section |lVF . So we only state the result: 



p-5G^(x,p) 



7r5(x^) sin [2x ■ p) — T'— cos (2a; • p) 



Now we define el projector onto the spRce perpendiculcir to the pa-rticle's momentum, — 9^i^ — Pfj-Pi' 

/p^. This 

allows us to change variables to a;_L^ — g_L^yX^ and t — x ■ pj \J\p'^\ sgn(p^). In terms of these variables, we find 
A^ = —p^x\ and the equation of motion becomes 



drG^ir, XA^,p) 



sgn(p2)yp 



'K5{\k^\T^ - A^) sin (2v/|fc2|T) - V 



cos (2v/|F|r) 

|fc2|T2-A2 



''The constraint -like equation could also be used, but is easier to derive using the transport-like equation. 
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So, instead of doing the Wigner transform directly, we only have to solve this ordinary differential equation. 

We find the solution by integrating this differential equation. The delta function integral is simple and the principle 
value integral can be done by contour integration. We find 



1 



2, sin(2VA2) 



e'(A' 



sgn(p2)0(_r)- 



2 



i9{p^) ^ e{x')) + sgn(p2)0(-r) 



We must now divine the boundary condition at t — > oo. 

To find the boundary condition, we actu ally have to go back to the Wigner transform of the propagator starting 
from momentum-space version of equation ( |GlD . We again change variable from a; to r and x± . We also change from 
p' to p'j^ = g±fj,i/p''^ and p-p' ^ sgn{p'^)^y\p^k. With this, we perform the k contour integral. The integral is straight 
forward, but tedious. However, when we take the limit as t — > oo, the result simplifies dramatically: 



G"(oo,a;_L,p) = -ri^ / d^P± cos (2a;_L ■p±)S{p^ +pl). 

\ \V \ J 



\v 

The delta function integral is trivial and the last pair of integrals requires integral tables, but in the end we find: 



Plugging this into the solution of our differential equation, we find equation ( p4bD . This result can be checked by 
performing a series of contour integrals in momentum or coordinate space. 



2. Massive Scalar Feynman Propagator 



Remlcr [ p^ has found the Wigner transform of the massive Feynman propagator. This transform is difficult but, 
when one uses the approximation p' « [p-p'Y /p^ , where p is the average momentum and p' is the relative momentum, 
the integrals become simple contour integrals. We state Rcmler's result here: 

drS^ix 



|-\/— p2 cos (2t-\/— p2) + msin (2|r|-\/— 



G%x,p) = { 



dT !-j= 

2m-y/p2 



sin (2r(-^p2_j^jj 
(-y/p^— m) 



for p"^ < 



for p2 > 



(G5) 



Note that, because of the approximation made, this propagator is oversmoothed in the direction transverse to the 
particle's momentum and we expect these propagators to be accurate only for length scales much larger than the 
size of the smoothing. Since the resulting propagators vary on length scales of order 1/m, we should only use these 
propagators for momenta with p m. Note also that the sine and exponential functions in the two terms in ( |G5| ) have 
the property that they become proportional to S(jP — m^) as r — > oo. Thus, this propagator reduces to the classical 
propagator Finally, we note that the 5-functions constrain the particle to more along its classical trajectory, even 
though its four-momentum (and hence its four-velocity) is being modulated by the sine and exponential functions. 
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TABLES 



TABLE I. The vertex Feynman rules for scalar and spinor QED. 



3 point 
photon-scalar 
vertex 




eZ *dy,= eZ(d^ - a^) 


4 point 
photon-scalar 
vertex 






fermion-photon 
vertex 


T 





TABLE II. The contour scalar, photon, and electron propagators. 



scalar line 




G{X1,X2) 


photon line 




Dft„{xi,X2) = 4:7Tgi^^G{xi,X2) 


fermion line 


^ 


Sal3{xi,X2) = -((/+ m)al3G{xi,X2) 
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FIGURE CAPTIONS 



FIG. 1. (a) Cut diagram for current A to exchange a photon with current B. (b) Cut diagram for current B to absorb a free photon. 
In both figures, the photon/current B interaction is unspecified and is represented with a blob. 

FIG. 2. Both figures arc plots of the dimcnsionlcss function A corresponding to tlic effective plioton distribution of a point cliargc with 
3-velocity v = (v]^,Ot) with = 0.9c. The photons in these slices of the phase-space distribution have = (me, rrie/i'L, Oy) (left) and 
q^l = {me, rUe/vL, 0.56 MeV/c, 0) (right). In both plots, only the negative and zero contours are labeled. The positive contours increase 
in increments of 0.25. 

FIG. 3. (a) Cut diagram for creating an electron-positron pair by photon splitting. The electron interacts with the probe particle, B. 
The square vertex represents the photon source, (b) Cut diagram for a free electron interacting with the probe particle. 

FIG. 4. On the left: the electron source for electrons with momentum = (2.0, 2.05, Ot) MeV/c electrons. In this figure, only the zero 
contours are labeled. The positive contours are (in arbitrary units) 1.0, 2.5, 5.0, 7.5 and 10.0. On the right: the virtual photon distributions 
corresponding to one of dominant contributions to electron source. These photons have {q+ij.) = (0.956,1.063,0.045,0.045) MeV/c. The 
other root has similar momentum and a similar distribution. In this figure, only the negative and zero contours are labeled. The positive 
contours increase in increments of 0.25 (in arbitrary units). 

FIG. 5. Coordinate space distribution of = (2.0, 2.05, 0) MeV/c electrons. Only the negative and zero contours are labeled. The 
positive contours are in increments of 1.0 (in arbitrary units). The sign of the contours in each region are denoted by ± signs. 

FIG. 6. Cut diagram for lepton pair production from a two photon interaction. R is the space-time point of the center of the collision 
region. 

FIG. 7. The diagrams that contribute, at lowest order, to the 77 ee effective vertex. 

FIG. 8. The ellipses represent the edge of the photon distributions, each with four-momentum q = {me,me/vi^,OT)- The shaded region 
is the geometrical overlap of the photon distributions and sets the size of the ee production region. The arrows point in the direction of 

the photons' source's 3— momentum. 

FIG. 9. The contour in the complex time plane used in the evaluation of operator expectation values. The upper branch corresponds 
to causal ordering and the lower branch to anticausal ordering. The arrows denote the contour ordering enforced by the T operator. 

FIG. 10. The Dyson-Schwinger equations of the propagators. Double lines represent the dressed Green's functions and single lines 
represent the non-interacting Green's functions. The particle self-energies are the large square vertices. 

FIG. 11. The scalar and electron self energies and the photon polarization tensor. Bare vertices are represented by dots and dressed 
vertices by blobs. The self-energies and the polarization tensor are all represented by large square vertices. 

FIG. 12. Cut diagram for probing the particle densities in the Generalized Fluctuation-Dissipation Theorem. Time flows downward 
and, since the the probe interaction is in the future, we leave the nature of the probe unspecified. 

FIG. 13. Cut diagram for the time-ordered (noncquilibrium) photon density. Time flows downward and, since the probe interaction is 
in the future, we leave the nature of the probe unspecified. 
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FIG. 14. Cut diagram for the time-ordered (nonequilibrium) electron density. Time flows downward and, since the the probe interaction 
is in the future, we leave the nature of the probe unspecified. 

FIG. 15. Cut diagram for probing the n*'' generation of partons in a typical cascade. Time flows downward in this diagram and the 
probe, being somewhere in the future, is left unspecified. 



FIG. 16. Typical rung of the LLA ladder. 

FIG. 17. Plot of the dimensionless function A. corresponding to the Wigner transform of the Coulomb field of a static point charge. 
The photons in this plot have = (0, 0.788, Ot) MeV/c. The longitudinal axis is defined by the photon three-momentum. 

FIG. 18. On the left: coordinate space distribution of spacc-likc (p^ = (0.05,0.008,0.06,0) McV/c) electrons. The ±'s indicate the 
sign of the function in a particular region. The contours go in steps of 25 (in arbitrary units). On the right; the photon distribution 
for photons with q+f^ = (-6.63,-6.65,0.0648,0.00478) MeV/c. The other root has q-^^ = (-6.33,-6.36,0.0648,0.00478) MeV/c and its 
distribution is similar. Again, the ±'s indicate the sign of the function in a particular region. Here, the contours go in steps of 0.25 (in 
arbitrary units). 

FIG. 19. Coordinate space distribution of time-like (p^ = (0.05,0.005,0.04,0) MeV/c) electrons. The ±'s indicate the sign of the 

function in a particular region. The contours arc (in arbitrary units) 150, 10, 5, 1, 0.2, 0, -1, -5, -10, and -150. One of the roots of the 
underlying photon distribution is shown at the right. These photons have momentum = (—8.39,-8.42,0.182,0.142) MeV/c. Here, the 
±'s indicate the sign of the function in a particular region and the contours go in steps of 0.25 (in arbitrary units). 
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